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I.  INTRODUCTION 


Most  current  vulnerability/lethality  programs  represent  a  penetrator's  path 
through  a  target  by  a  shotline.  This  shotline  is  constructed,  at  least  conceptually,  by 
passing  a  ray  through  a  geometric  description  of  the  target  and  recording  each  target 
component  that  the  ray  intersects  and  at  what  distance  from  its  origin  the  ray  enters 
and  exits  the  component.  The  analyst  will  typically  select  azimuth  and  elevation 
angles  from  which  to  view  the  target,  and  then  superpose  over  the  view  of  the  target 
(hereafter  referred  to  as  V)  a  rectilinear  grid,  as  shown  in  figure  1.  From  a  selected 
location  with  each  cell  in  this  grid  a  ray  will  be  sent  through  the  target  to  generate  a 
shotline.  The  aggregate  of  these  shotlines  is  then  used  in  determining  the  entire 
target’s  vulnerability. 

One  drawback  of  the  shotline  approach  is  that,  whereas  many  man-months 
might  be  spent  on  a  detailed  description  of  the  geometry  of  the  target,  the  projectile 
can  only  be  represented  by  a  straight-line  trajectory,  an  abstraction  with  zero  cross- 
sectional  area.  Of  course,  if  the  projectile’s  geometry  and  fuze  design  are  accommo¬ 
dating,  and  in  the  absence  of  yaw  and  other  complications,  this  limitation  is  not  too 
serious  for  planar  targets  of  infinite  extent.  But  for  some  projectiles  —  either  with 
piezoelectric  contact  fuzes  on  probes  or  with  inductance-  or  capacitance-type  influ¬ 
ence  fuzes  —  the  single-ray  approach  is  insufficient.  Even  such  common  target 
characteristics  as  discontinuities  and  high  obliquities  cause  problems  for  zero-width 
modeling  of  projectiles. 

This  report  presents  a  modification  of  the  standard  shotline  method.  A  com¬ 
puter  routine  called  SHOTCYL  provides  a  means  of  modeling  the  three-dimensional 
geometry  of  projectile/target  interaction.  The  projectile  is  represented  by  a  main,  or 
central,  ray  together  with  one  or  more  rings  of  planetary  rays  disposed  parallel  to 
and  at  specified  distances  from  the  main  ray. 

This  representation  has  been  used  in  analyzing  the  behavior  of  several  HEAT 
rounds.  In  such  a  warhead,  the  charge  is  mounted  some  distance  behind  the  fuze. 
This  built-in  standoff  is  intended  to  allow  the  penetrating  jet  to  form  before  it 
strikes  the  target.  A  jet’s  penetration  into  armor  is  sensitive  to  standoff,  so  any 
abnormal  impact  (such  as  that  illustrated  in  figure  2)  might  significantly  degrade  a 
round’s  performance. 

Given  a  bundle  of  rays  for  each  grid  cell,  a  program  called  MISFIR  computes 
the  effective  standoff  for  each  cell,  using  the  built-in  standoff  and  the  geometrical 
details  of  the  projectile’s  impact  on  the  target.  This  provides  for  more  realistic 
treatment  of  several  projectile  designs  and  of  such  target-surface  properties  as  obli¬ 
quity  and  edges. 
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Figure  1.— View  of  a  Target  with  Grid  Superposed. 
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Fiovus  2.— Edge  Hit  of  a  Shaped  Charge  Projectile 
on  an  Armored  Target 


Also  discussed  here  is  the  employment  of  the  new  method  in  an  analysis  of  the 
effectiveness  of  a  hypothetical  two-fuzed  missile.  The  missile  has  a  primary  fuze 
that  is  mounted  on  a  long  probe  and  a  secondary  fuze  further  back  on  the  shoulder 
of  the  missile,  as  shown  in  figure  3.  The  modified  shotline  method  was  used  to 
answer  two  questions: 

1)  For  a  given  projectile,  target,  view,  and  aimpoint,  what  is  the  probabil¬ 
ity  that  the  secondary  fuze  will  strike  the  target  before  the  primary  fuze? 

2)  Under  the  same  conditions,  what  is  the  probability  that  the  secondary 
fuze  will  strike  the  target  after  the  primary  fuze,  but  during  the  primary’s 
built-in  delay? 


II.  PROJECTILE  MODELING  TECHNIQUE 

The  key  input  for  the  enhanced  method  is  a  description  of  the  projectile. 
Depending  on  the  nature  of  the  projectile,  this  three-dimensional  model  will 
represent  the  projectile’s  fuze,  its  nose,  its  influence  envelope,  or  whatever  the 
relevant  volume.  The  following  paragraphs  provide  a  conceptual  explanation  of  the 
model. 

Consider  an  arbitrary  solid  S  traveling  along  a  straight-line  trajectory,  T.  The 
volume  it  sweeps  out  is  a  cylindroid.  We  add  a  frame  of  reference  whose  origin,  O, 
is  at  the  front  tip  of  S.  Now  let  us  make  the  assumption  that  the  solid  is  radially 
symmetric  about  T,  and  that  its  maximum  circumference,  C,  occurs  at  some  distance 
D  along  T  (see  figure  4).  A  snapshot  (ignoring  the  early  part  of  the  trajectory) 
reveals  that  the  shape  of  the  volume  swept  out  is  that  of  a  right  circular  cylinder 
capped  with  the  patch,  P,  of  S' s  surface  that  is  bounded  by  C  and  contains  O. 

Any  collision  between  S  and  a  stationary  object  will  take  place  at  a  point  on  P. 
This  patch  is  the  crucial  part  of  our  model,  and  can  be  approximated  by  a  set  of  cir¬ 
cles  that  are  centered  on  7,  each  circle  being  specified  by  its  radius  and  its  stepback 
—  the  distance  from  O  to  the  circle’s  center.  As  the  solid  travels,  the  circles  sweep 
out  concentric  cylinders.  These  cylinders,  as  they  exist  at  any  specified  instant,  can 
be  represented  by  rings  of  rays  originating  on  the  circles  and  extending  backwards. 

In  principle  this  is  how  a  projectile  is  modeled.  The  program  determines  at 
what  point  each  ray  intersects  the  target  and,  considering  the  contours  of  the 
projectile’s  leading  surface  (as  rendered  by  the  model),  decides  where  and  how  the 
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Figure  3. — Hypothetical  Two-Fuzed.  Missile. 


Figure  4. — Flight  of  a  Radially  Symmetric  Solid,  S, 
Modeled  as  a  Bundle  of  Rays. 
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projectile  first  touches  the  target.  This,  together  with  the  projectile’s  built-in  stand¬ 
off,  provides  enough  information  to  compute  the  effective  standoff  for  the  hit  being 
considered. 


-TO.  SYNOPSIS  OF  ALGORITHM 

The  three-dimensional  model  of  projectile/target  interaction  is  centered  upon 
the  program  MISFIR,  written  in  CDC  Fortran  5.  MISFIR  is  built  on  the  formalisms 
of  the  GIFT  (Geometric  Information  /or  Targets)*  program,  which  processes  targets 
described  in  terms  of  combinatorial  geometry,  or  COMGEOM.  The  MISFIR  pack¬ 
age  currently  consists  of  a  ray-tracing  subroutine  added  to  GIFT  (viz.  SHOTCYL); 
MISFIR  itself,  together  with  its  subprograms;  and  an  application  program,  called 
FUZES,  which  uses  MISFIR’s  results  to  solve  a  typical  problem  in  vulnerability 
analysis.  A  generalization  of  the  familiar  GIFT  shotline  provides  the  means  by 
which  MISFIR  represents  projectile  geometry:  The  single  ray  is  augmented  by  a 
bundle  of  rays,  the  combination  being  called  a  shotcylinder.  The  planetary  rays  are 
grouped  around  the  central  ray  in  orbits  that,  depending  on  their  radii,  can  extend 
into  adjacent  grid  cells. 

The  first  step  is  to  create  shotcylinders  for  each  cell  in  V,  the  current  view  of 
the  target.  This  is  done  by  the  GIFT  subroutine  called  SHOTCYL.  Next  the 
shotcylinders  are  used  to  compute  the  projectile’s  actual  standoff  in  each  cell.  MIS¬ 
FIR  then  creates  histograms  of  the  standoff  and  a  silhouette  image  of  the  target.  In 
the  application  described  here,  FUZES  is  then  run  to  analyze  the  effectiveness  of  the 
hypothetical  warhead.  Below  are  brief  descriptions  of  each  step  in  the  process. 

The  user’s  input  to  SHOTCYL  includes  the  angle  of  attack,  expressed  in  terms 
of  azimuth  and  elevation,  and  such  ray-bundle  building  parameters  as  the  number  of 
orbits  [layers]  per  bundle,  the  orbits’  radii,  and  the  number  of  planets  [rays]  in  each 
orbit.  The  frame  of  reference  for  GIFT’S  calculations  is  defined  by  the  grid  plane, 
the  plane  that  contains  the  origin  of  the  target  description’s  coordinate  system  and  is 
normal  to  the  line  of  sight,  or  attack  angle.  The  grid  plane  is  partitioned  into  rectil¬ 
inear  cells,  and  the  target’s  projection  onto  it  determines  the  cells  for  which  bundles 
must  be  created. 


*  Kuehl,  G.G.,  Bain,  L.W.,  Jr.,  and  Reisinger,  M.J.  The  GIFT  Cod e  User  Manual. 
BRL  R1802,  July  1975,  and  ARBRL-TR-02189,  September  1979. 
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Also  input  to  SHOTCYL  is  the  manner  of  selecting  the  central  ray’s  location 
within  the  grid  cell.  SHOTCYL  can  use  an  impact-point  file  to  treat  a  target  that 
has  already  been  shotlined*.  Alternatively,  SHOTCYL  can  pass  its  central  rays 
through  either  the  center  of,  or  a  point  chosen  at  random  in,  each  grid  cell.  For 
each  grid  cell,  SHOTCYL  outputs  the  following  information: 


variable  names 

contents 

HCEN,  VCEN 

The  grid-plane  coordinates  of  the  center  of  the  cell 

IHIV 

A  code  for  the  location  of  the  central  ray 
within  the  cell 

H,  V 

The  grid-plane  coordinates  of  the  central  ray 

RAYLEN 

The  distance  along  the  central  ray  from 
the  grid  plane  back  to  the  first  item  met. 

Thus,  this  number  is  negative  for  cells  in 
which  the  central  ray  does  not  encounter  the 
target  until  it  has  passed  through  the  grid  plane. 

IT M HIT 

The  item  number  of  the  first  item  met  by  the 
central  ray,  or  -1  if  it  never  hits  the  target. 

And  for  each  planetary  ray: 

PLNETH,  PLNETV 

The  grid -plane  coordinates  of  the  ray 

RAYLEN 

The  distance  along  the  ray  from  the  grid  plane 
back  to  the  first  item  met.  Thus,  this  number 
is  negative  for  planetary  rays  which  do  not 
encounter  the  target  until  they  have  passed 
through  the  grid  plane. 

TTMHTT 

The  item  number  of  the  first  item  met  by 
the  ray,  or  -1  if  it  never  hits  the  target. 

High  resolution  in  shotlining  the  typical  target  is  not  always  computationally 
feasible.  The  relatively  large  grid  cells  resultant  from  this  constraint  render  many 
applications  programs  very  sensitive  to  ray  location  within  a  cell.  In  order  to 
control  for  positional  variation,  Robert  Wilson  has  developed  a  technique  of 
recording  cell  impact  points  in  a  file  parallel  to  the  shotline  file.  When  a  target 
must  be  shotlined  repeatedly,  as  in  a  parametric  analysis,  the  ray  coordinates  for 
each  cell  can  be  reused,  thus  maintaining  comparability. 
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SHOTCYL  requires  one  subprogram  that  is  not  found  in  the  GIFT  or  system 
libraries.  It  is  a  routine  called  SEEKVEW  that  scans  through  an  impact-point  file  to 
find  the  data  for  a  specified  target  aspect. 

MISFIR  is  the  centerpiece  of  the  package.  It  produces  its  results  —  information 
about  projectile  impact  for  each  cell  on  the  target  —  in  several  forms.  The  subpro¬ 
grams  called  by  MISFIR  are  CPA,  PHIT,  HISTOG,  and  SILOET. 

CPA  selects  as  the  aimpoint  the  target’s  center  of  presented  area.  This  is 
equivalent  to  the  centroid  of  V,  and  can  be  calculated  as: 


where  N  is  the  number  of  cells,  c,  in  V.  This  aimpoint  is  then  offset  by  ( Xdc ,  Ydc)  to 
allow  for  dispersion  corrections  provided  by  the  user. 

Once  the  aimpoint  has  been  computed,  MISFIR  determines  which  ray  in  each 
cell  strikes  the  target  first  and  calculates  the  size  of  the  gap  between  projectile  tip 
and  target  at  the  instant  of  impact  (see  figure  5).  The  basic  formula  for  this  gap  is: 

GAP  =  max  (0  ,  max  (DPLNETp  -  DMA1N  -  CSTPBKp  ~  CSTPBK0))  (2) 

p  in  B 

where  p  is  a  planetary  ray  in  B ,  the  current  cell’s  bundle,  DPLNETp  is  the  distance 
along  p  from  the  grid  plane  back  to  the  first  planet  met,  DMAIN  is  the  distance 
along  the  central  ray  from  the  grid  plane  back  to  the  first  item  met,  CSTPBKp  is  the 
stepback  of  p’s  cylinder,  and  CSTPBK0  is  the  length  of  a  needle  probe  of  negligible 
diameter. 

For  cells  in  which  the  primary  fuze  strikes  the  target  before  the  secondary  fuze, 
MISFIR  also  determines  whether  the  secondary  fuze  will  strike  during  the  primary 
fuze’s  delay.  Then  MISFIR  determines  the  warhead’s  actual  standoff  for  that  cell  — 
the  projectile’s  built-in  standoff  plus  GAP  —  and  the  probability  of  the  cell’s  being 
hit  (Pj/).  The  cell’s  P#  is  computed  by  the  function  PHIT,  which  assumes  that 
delivery  errors  are  normally  distributed. 

HISTOG  computes  and  displays  two  histograms.  To  do  this,  it  partitions  into 
bins  the  range  of  standoffs  that  were  encountered  in  V.  The  first  histogram  gives 
the  number  of  cells  whose  standoffs  lie  within  each  bin.  The  second  displays  similar 
data  obtained  by  choosing  an  aimpoint  on  the  target  and  using  the  weapon’s 
delivery-accuracy  characteristics  to  weight  the  standoff  obtained  in  each  cell  by  that 
cell’s  P# .  Thus,  the  histograms  convey  the  relative  frequencies  of  various  values  of 
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DPLNETi 
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Figure  5. — Geometry  of 


standoff,  either  assuming  a  random  hit  anywhere  on  the  target  or  assuming  a  shot 
aimed  at  its  centroid.  Both  histograms  also  include  a  bin  for  edge  hits,  cells  where 
at  least  one  planetary  ray  intersects  the  target,  but  the  main  ray  misses  completely. 
SILOET*  displays  an  image  of  the  non-empty  portion  of  the  grid  plane,  where  each 
cell  is  represented  by  standoff  to  two  decimal  places,  and  in  units  of  tenths  of  a 
charge  diameter.  MISF1R  also  prints  all  its  results  to  a  file  for  use  by  other  pro¬ 
grams.  The  information  produced  for  each  cell  is: 


variable  names 

contents 

HCEN,  VCEN 

The  grid-plane  coordinates  of  the  cell. 

DIST 

The  effective  standoff  for  the  cell. 

CYL1 

The  number  of  the  cylinder  containing  the  ray  that 
struck  the  target  first. 

FZ2HTT 

A  Boolean  flag  indicating  whether  both  the  primary 
fuze  struck  the  target  first  and  the  secondary 
fuze  struck  during  the  primary  fuze’s  delay. 

CELLPH 

The  cell’s  P#. 

FUZES  provides  information  about  the  performance  of  the  hypothetical  missile 
by  estimating  the  conditional  probabilities  of  the  two  events 

E\.  The  secondary  fuze  strikes  the  target  before  the  primary  fuze  does. 

and 

£2:  The  secondary  fuze  strikes  the  target  after  the  primary  fuze  does, 
but  during  the  primary  fuze’s  built-in  delay. 

both  given  a  uniform-random  hit  on  the  target,  and  given  that  a  shot  was  fired  at  the 
specified  aim  point.  Computing  the  probabilities  of  Ex  required  only  one  modifica¬ 
tion  to  MISFIR:  the  value  of  an  existing  variable  (viz.  CYZ.1)  would  have  to  be 
added  to  the  output.  £2  induced  a  slightly  larger  task.  Several  lines  of  code  had  to 
be  added  to  MISFIR  to  handle  the  computation  and  the  printing  of  the  value  of 
FZ2HIT. 


Based  on  H.  Ege’s  SILPK  program 
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The  probability  of  E\  given  a  hit,  P(£j  |  H),  is  simply  the  ratio  of  the  number 
of  cells  in  which  the  secondary  fuze  struck  first  (i.e.  the  number  of  cells  for  which 
CYL\  exceeded  some  user-specified  threshold)  to  the  the  number  of  cells  in  V.  The 
probability  of  E\  given  a  shot,  P(£i  (5),  is  the  sum  of  the  Pff’s  of  all  the  cells  in 
which  the  secondary  fuze  struck  first. 

FUZES  computes  the  probability  of  £2  given  a  hit  by  using  the  formula 

P(£2  Iff)  =  2  Fzltirr-  (3) 

c  in  V 


Similarly,  P(£2  |  S)  is  the  sum  of  the  Pw’s  of  all  the  cells  for  which  FZ2HTT  =  1. 

The  source  code  for  SHOTCYL,  MISFIR,  FUZES,  and  their  subprograms  can 
be  found  in  appendix  A.  Appendix  B  provides  a  sample  run  of  the  package.  The 
histograms  on  pages  60  and  61  indicate  that,  whether  the  hits  were  distributed  uni¬ 
formly  or  normally,  the  preponderance  of  hits  occurred  at  built-in  standoff.  They 
also  show  that,  although  there  was  about  a  13%  probability  of  achieving  an  edge  hit 
given  a  random  hit,  when  delivery  accuracy  was  considered,  the  probability  was 
reduced  by  an  order  of  magnitude.  The  silhouette  on  page  62  illustrates  these  points 
for  each  cell.  For  the  numbers  in  the  sample  run,  the  very  common  silhouette  cell 
value  of  25  represents  2.5  charge  diameters,  or  250  mm,  which  is  the  projectile’s 
built-in  standoff.  FUZES’  output,  page  63,  indicates  that  P(E{  |  H)  ~  22%,  and 
P(£j  |  S)  ~  12%,  while  P(£2  |  H )  is  about  8%,  and  P(£2  |  S )  is  about  5%. 


IV.  CONCLUSIONS 

Initial  runs  of  the  SHOTCYL-MISFIR-FUZES  package  pitting  the  hypothetical 
missile  against  one  important  target  suggest  that  the  likelihood,  and  therefore  the 
importance,  of  secondary-fuze  detonation  can  be  quite  significant.  More  thorough 
testing  —  considering  several  projectiles  and  targets  —  could  perhaps  settle  the  issue 
more  conclusively.  In  any  event,  the  MISFIR-FUZES  model  should  answer  the 
question  easily  and  cheaply.  It  ought  therefore  to  be  useful  for  designers  of  war¬ 
heads  and  armored  systems,  and  for  those  of  us  who  assess  the  vulnerability  of 
existing  systems.  BRL  will  be  using  MISFIR  routinely  in  P*  analyses. 
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Additions  and  modifications  to  the  program  will  no  doubt  be  required  in  time. 
Some  changes  that  are  already  under  consideration  are: 

1)  To  allow  the  number  of  planetary  rays  per  orbit  to  vary  within  the 
bundle. 

2)  To  allow  the  option  that  planetary  rays  be  positioned  at  random 
around  the  orbits.  MISFIR  currently  spaces  planetary  rays  evenly  around 
the  orbit. 

3)  To  change  the  format  of  the  histogram  output  for  clarity  and  more 
precision. 
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APPENDIX  A 
Source  Listings 
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oooooooooonoonoonooonnoonoor^ononor^ooriooooooooooooooooooriri 


SUBROUTINE  BRANDX 


A.X.A.  SHDTCYL 


COMPUTES  FIRST  INTERSECT  FOR  CYLINDERCS)  AROUND  A  MAIN  GRIO  RAY 


DERIVED  FR  DM  GART  KUEHL'S  'GRID'  SUBROUTINE. 
SIGNIFICANTLY  EMBELLISHED  BY  PAUL  TANENBAUM 


CONTPDL  CARD  (315) 
1-5  NMASPS 

5-10  NDPRNT 
11-15  RAYLOC 


NUMBER  OF  VIEWS  FOR  CYLINDER 
(CURRENTLY  MUST  BE  1) 

.NE.  0  *•>  SUPPRESS  PRINT  OF  UNIT  1  DATA 
LOCATION  OF  RAT  WITHIN  GRID  CELL 
(.EQ.  0  ==>  RANDOM  POINT  IN  CELL) 

(.SO.  1  **>  READ  FROM  RAYFILE) 

C.GE.  2  =■>  CENTER  OF  CELL) 


PROJECTILE  card  CA23, 
1-20  RJNNAM 
21-25  NH-CrLS 

25- 3 G  NPlnTS 


31-37  CD 
3v-4<»  STDOFF 
45-S2  VEL 
53-50  DELAY 
61-65  *Z1MXC 


215.  2F7.2,  2F8.3.  15) 

-  CHAR  STRING  TD  IDENTIFY  RUN 

-  NUMBER  OF  CYLINDERS 

(1  <-  NWCYLS  <  =  10.  DEFAULT  «■  2) 

-  NUMBER  OF  RAYS  PER  CYLINDER 

(1  <=  NPLNTS  <■  50,  DEFAULT  »  B) 

-  CHARGE  DIAMETER  CNN 3 

-  BUILT-IN  STANDOFF  CNN) 

-  VELOCITY  CM/SI 

-  BUILT-IN  OELAT  OF  PRIMARY  FU2E  C  SEC  3 

-  NUMBER  OF  LARGEST  CYLINDER  IN  PRIMARY  FUZE 


PROJECTILE- SKIN  CARDS  (3F10.2)  CCNE  CARO  PER  CYLINDER  3 
1-10  CRAS(N)  -  RADIUS  OF  NTH  CYLINDER  CMH3 

11-23  CSTP-5CN)  -  STEP6ACK  DF  NTH  CYLINDER  CNM3 

21-30  CSTPi((D)  -  LENGTH  OF  Z £ R 0-C ROS S-SE C TI DNAL- A R E A 

PROBE  CNH3 

CLAST  PROJECTILE-SKIN  CARD  DNLt) 


CIS’ERSIDN  CARD  (4F3.3) 
1-3  SI GX 
3-16  SIS Y 
1 7- 2  A  XDC 
25-32  Y DC 


X-DI $°ER5IQN  CMN3 
Y-DISPERSION  CNN) 

X  DISPERSION  CORRECTION 
T  DISPERSION  CORRECTION 


CNN  3 
CHM3 


ASPECT  CARDS  (3F1D.O, 
1-10  AZIM 
11-20  eLEV 
21-30  CELSIZ 

31-35  IDLKP 
36-40  1LIM1T 
*1-45  ISPOT 


315)  CONE  CARD  PER  VIEW) 

-  AZIMUTH  ANGLE  CDEGREES3 

-  ELEVATION  ANGLE  C0EGREES3 

-  SUE  DF  CELLS  IN  GRIO  PLANE  CMH3 
(DEFAULT  ■  *  IN) 

-  DELETE  SELECTED  REGIONS 

-  LIMIT  GRID  PLANE 

-  NUMBER  OF  SELECTED  SHOTLINES 


LIMIT  GRID  PLANE  CARD  (AF10.0)  CIF  ILIMIT.NE.03 

1-10  HMIN  -  HORIZONTAL  MINIMUM  OF  GRID  PLANE 

11-20  HM AX  -  HORIZONTAL  MAXIMUM  OF  GRID  PLANE 

21-30  VMIN  -  VERTICAL  MINIMUM  OF  GRID  PLANE 

31-40  V M A X  -  VERTICAL  MAXIMUM  OF  GRID  PLANE 
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C  GRID  PLANE  COORDINATES  CARDS  (3-10.0)  CIFF  ISPCT  .GT.  03 

C  1-10  CMC1)  -  HORIZONTAL  COORDINATE  OF  POINT  l 

C  11-20  C V C 1 3  -  VERTICAL  COORDINATE  OF  POINT  1 

C 

C  71-80  C  V  ( 4 )  -  VERTICAL  COORDINATE  OF  POINT  A 

C 

Character  tio  day 
CHARACTER  *60  ITITLE 
CHARACTER  *20  SUNNAN 

INTEGER  ASPECT,  CTL,  FIIHXC,  I,  IOLXP,  IDUM,  IHIV,  IHDRZ,  ILIHIT, 

1  INPCEL,  IRANH,  I RANV,  IREGDN,  ISEED,  ISEED1,  ISOLID. 

2  ISPCT,  ISL'RF,  ITHHIT.  IVERT,  J,  L,  LENRPP,  L I  RFC  t  NCELL, 

3  NHORZ,  NHASPS,  NMCYLS ,  NHHITS,  NOPRNT,  NPLNTS.  NVSRT, 

a  planet,  ratldc 

LOGICAL  HK..3AYS,  PRT ALL 

REAL  A3  SC A ,  ASSCE,  AESSA,  ASSSE,  ANGINC,  ASTER,  AZIH,  A  Z I  HR , 

1  6ACK.,  CAZIH,  CD,  CELEV,  CELSIZ,  CM,  CRAD,  CSTP9K,  CV, 

2  CEG2RAD,  DELAY,  DEPTH,  0  F  I  3  S  T  ,  ELEV,  ELEVR,  ENGTH,  EOV, 

3  ETIHE,  GCLIH,  H,  HC,  HCEN,  MMAX,  HH  IN ,  HDRZ,  HMCVT,  PLNSTH, 

a  PLNETV,  RAYLEN,  RFA,  RFCS,  RFE,  SANGLE,  SAZIH,  SELEV,  SIGX, 

5  SIGY,  SJIGlE,  STCDFF,  stihe,  T,  TCEN,  TGTCVT,  TLEN,  THAX, 

6  THIN,  V,  VC,  VCEN,  VEL,  VERT,  VHAX,  VHIN,  VTIHE,  W3,  We,  *6, 

7  XuC ,  YOC 
LEVEL  2,// 

CCHHDN  /ENSP?  /  LEN  S  P  0 

CCHHCN  /GECH  /  GDUH(a),  LIRFO 

CCMHDN  /DNEINT/  DF 13  ST ,  ISOLID,  ISURF,  IREGDN 

COMMON  / 3  AT  PA  R/  X2(3),  trSC3) 

COHHCN  /SEED  /  ISEED 
COMMON  /TGTCVT/  TGTCVT 
COMMON  /TITLE  /  ITITLE 

DlNENSICN  CH(4),  C3ADC10),  CSTP9K.C 0:10),  CVCA), 

1  I TMHIT  (0 :  10 ,  0:50),  RATLEN(0:lO,  0:50),  PLNETHOO,  50), 

2  PLNETVOO,  50),  T*IN(3),  TMA  x  ( 3 )  •  TCEn<3).  TLENC3), 

3  RO(3) 

COMMON  ASTERC5000) 

OATA  OEG2RAD  /. 0 1745329251 9943/ 

DATA  SJIGLE  /.0001/ 

C  -  SEGIN  EXECUTION  - 

HHCVT  «  25. A  /  TGTCVT 

C  -  TGTCVT  -  1  IF  TARGET  DESCRIPTION  IS  IN  INCHES, 

C  »  25.4  IF  TARGET  DESCRIPTION  IS  IN  NM. 

C  THEREFORE,  DIVIDING  A  LENGTH  IN  HM  3 Y  HHCVT  HILL  CONVERT 

C  THE  LENGTH  TO  THE  UNITS  OF  THE  TARGET  DESCRIPTION. 

C  CONVERSELY,  MULTIPLYING  A  LENGTH  IN  THE  DESCRIPTION'S 

C  -  UNITS  6T  HHCVT  WILL  CONVERT  THE  LENGTH  TO  NM. 

call  CLOCKS(STIHE)  ■ 

c -  RE  AO  CONTROL  CARD 

READ  (»,  5010,  ERR-160,  END-170)  NHASPS,  NOPRNT,  RATLOC 
C  — -  NHASPS  MUST  BE  .ED.  1...  IF  NOT,  ABORT 

IF  (NHASPS  .NE.  I)  THEN 
WRITE  (*,  6010) 

STOP 
END  IF 


26 


C -  DEFINE  PRT  ALL ,  AND  IF  IT'S  FALSE.  WARN  AS  MUCH 

PRTALL  •  NCPRNT  .EC.  0 
IFC.NCT.  PRTALL)  THEN 
WRITE  (P,  602  0 ) 

ENDIF 

C -  READ  PR  D JECTILE  DATA 

REAL  (P,  5020  .  ERR=160,  END-170)  RJNNAH,  NHCTLS.  NPLNTS.  CD. 

1  STDCFF,  VEL.  DEL  AT.  FI1MXC 

C -  ENSURE  THAT  FZ1PXC  IS  N  3  LARGER  THAN  NHCTLS 

IF  (FIIMXC  .GT.  NHCTLS )  THEN 
WRITE  C«.  6030)  FI1NXC,  NHCTLS 
STOP 
ENDIF 

C  -  ENSURE  THAT  1  <=  NHCTLS  <=  10 

IF  CNHCTLS  .LE.  0)  THEN 
NHCTLS  c  2 
ELSE 

NHCTLS  «  HINCNHCTLS,  10) 

ENDIF 

C  -  ENSURE  THAT  1  <=  NPLNTS  <=  50 

IF  CNPLNTS  .L£.  C)  THEN 
NPLNTS  *  6 
ELSE 

NPLNTS  *  HINCNPLNT5,  50) 

ENDIF 

C -  READ  CYLINDER  RADII  AND  STS  ®9 ACX  S. 

C  -  IF  RADII  ARE  NOT  POSITIVE,  IN  NON-DECREASING  ORDER.  ABORT. 

CO  10  CTL  *  1,  NHCTLS 

IF  CCTL  .Ei.  NHCTLS)  THEN 

READ  CP.  5030,  ERRsliO,  END  *170)  C8A3CCTL),  CSTP6KCCYL), 

1  CSTP3KC0) 

ELSE 

READ  CP,  5030,  SRR-150.  END-170)  CRADCCTL),  CSTPoKCCTL) 
ENCIF 

CRADCCTL)  *  CRADCCTL)  /  HHCVT 
CSTP3HCCTL)  =  CSTPiKCCTL)  /  MNCVT 
10  CDNTINJE 

CSTP3KC0)  «  CSTP5KC0)  /  HHCVT 
DO  20  CTL  ■  1.  NHCTLS  -  1 
IF  CCCRACCCTL)  .EC.  0)  .OR. 

1  (CRADCCTL)  .GT,  CRADCCTL  ♦  1)))  THEN 

WRITE  CP,  6040) 

STOP 

ENDIF 

20  CONTINUE 

C  -  DEFINE  HKRAYS 

HXRATS  3  CRATLCC  • N£,  1) 

C  — —  READ  DISPERSION  DATA 

READ  (P,  5040,  ERP-160,  EN0-1T0)  SIGX,  SIGT,  XDC,  T DC 
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CALL  DATECDAY) 
re-ino  i 

IF  C  PR  T  AL  L )  THEN 

WRITE  (*,  6050)  NMASPS,  OAT,  ITITLE 

WRITE  CP,  6060)  PSTALL,  RUNNAM,  CD.  STDOFF,  VEL,  DELAY,  FZ1MXC 
WRITE  (*,  6  0T  0  )  $  1 3  K  .  SIDY,  XDC,  YOC 
ENDIF 

WRITE  <1,  6050)  NMASPS,  DAY,  ITITlE 

WRITE  Cl,  6060)  PR T ALL ,  RUNNAM,  CO,  STDOFF,  VEL,  DELAY,  FZ1MXC 
WRITE  Cl,  607  0)  S I GX  ,  SIGY,  XDC,  Y DC 

C -  SEED  FOR  RANDOM  N'JMJfcR  GENERATOR 

ISEED  *  0 

C -  CYLINDER  RAY  TRAD- 

WRITE  Cp.  6060) 

DO  30  CYL  =  I,  NMCYLS 

WRIT;  (6,  6030)  CYL,  CRADCCYL)  *  MHCVT,  CYL,  CSTP3KCCYL)  *  MHCVT 
30  CONTINUE 

IF  CCSTP3XC0)  . j  T.  0)  TrEN 
WRITE  Cp.  6100)  CSTP3HC0) 
end:? 

WRITE  Cp,  6110)  NPLNTS 
IF  CPRTALL)  TnSN 

WRITE  C  9 ■  6120)  NMCYLS,  NPLNTS 
00  AO  CYL  *  1,  NMCYLS 

WRITE  Cp,  6130)  CST»3XCCTL) 

AO  CONTINUE 

•RITE  CP,  6130)  CST»3XC0) 

ENOIF 

WRITE  Cl,  6120)  NMCYLS,  NPLNTS 
DC  SO  CYL  »  1,  NMCYLS 

WRITE  Cl,  6130)  CSTP3KCCYL) 

50  CONTINUE 

WRITE  Cl,  6130)  CSTPSXCO) 

C  -  PROCESS  EACH  REQUESTED  VIEW 

DC  150  ASPECT  =  1,  N  M  A  S  p  S 
CALL  CLCCRSCVTIME) 

INPCEL  *  0 
ISEEDl  *  ISEED 

READ  CP.  5050,  EBpalSO,  END=170)  AZIM,  ELEV,  CeLSIZ,  IOLKP, 

1  ILIM1T,  ISPOT 

IF  CELSIZ  IS  NOT  GIVEN,  OR  IS  ZcRC,  SET  IT  TO  A  INCHES 
CCCNVERTcO  TO  Tnc  UNITS  OF  ThS  TARGET  DESCRIPTION). 

OTHERWISE,  SIMPLY  CONVERT  IT  7D  THE  UNITS  OF  THE  TARGET 
-  DESCRIPTION. 

IF  C C EL  SI 2  .IS.  0.)  THEN 
CELSIZ  «  4.  *  TGTCVT 
ELSE 

CELSIZ  -  CELSIZ  /  HMCVT 
ENDIF 
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IF  RATS  MUST  EE  MACE,  then  THEIR  LOCATION  HILL  Be  HARDWIRED  TO 
RANDOM  POINT  IN  GRID  CELL. 

JF  CMAS4T5)  THEN 
RAYLOC  «  0 
ELSE 

CALL  SEEK V£  W(51  *  'RAYFJLE*,  A  Z I M. ) 

R  E  A  j  (  51*  5360*  ERR*180,  END=19D)  RFA  ,  R't,  RFC5 
IF  (CELEV  •  N£  •  RF£)  .OR.  (CELSIZ  .NE.  RFCS))  THEN 
WRITE  («*  5140)  AZIN,  ELEV.  CELSIZ.  RFA*  RFE,  RFCS 
STOP 
£  NO  I F 
ENDIF 

IF  (IOLKP  .Nr.  0)  ThEN 
CALL  DELETE 
E.NOIF 

IF  (ISPCT  .GT.  0)  THEN 
RAIL  DC  *  -I 
ENDIF 

RETRIEVE  TARGET  MIN  AND  MAX 
L  =  L  ENR  -’P 
00531“  1*3 

TMIN(I)  «  ASTER(L) 

TMAX(I)  »  ASTERCL  ♦  1) 

TLEN(I)  *  TMAX(I)  -  TM.IN(I) 

TCEN(I)  «  .5  *  (TKAX(I)  ♦  TMISlCD) 

L  *  L  ♦  2 
C  ON T  I  Nil  E 

DIRECTION  COSINES 
A 1 1 M  A  I  All-  *  OEO’SAD 
ElEVR  =  ELEV  *  0EG2R4C 
SAZIM  *  SIN(AHMR) 

CAZIM  *  COKAZIM.R) 

SELEV  =  SIN(ELEVR) 

CELEV  =  COS (ELEVR) 

Wi(l)  »  -CELEV  *  CAZIN 
Us( 2 )  =  -CELEV  *  SAIIM 
Wa(3>  =  -SELEV 

COMPUTE  DIMENSIONS  OF  GRID  PLANE 
AES  S A  ■  AsS(SAZIM) 

ABSCA  *  AEKCAZIH) 

AbSSE  =  ASSCSELEV) 

ABSCE  »  A5SCCELEV) 

ENGTH  =  TLENC1)  *  ABSCA  ♦  TLENC2)  •  ABSSA 
HDRZ  ■  A5SCTLENC1)  *  ABSSA  ♦  TLENC2)  *  ABSCA) 

VERT  »  ABSCENGTH  *  ABSSE  ♦  HENO)  *  ASSCE  ) 

DEPTH  a  ABSCENGTH  *  ASSCE  ♦  TLENC3)  *  ABSSE) 

BACK  «  A1NT(0EPTH  *  .55  -  DDT(W&,  TCEN)  ♦  .5) 

FIND  COORDINATES  OF  CENTER  OF  GRID  PLANE 
T  ■  TCEN(l)  *  CA2IM  ♦  TCEN(Z)  *  SAZIM 
HC  •  -TCeN(l)  *  SAZIM  ♦  TCEN(2)  *  CAZIM 
VC  ■  -T  *  SELEV  ♦  TCENC3)  *  CELEV 
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IF  (ILIHIT  .£3.  0)  THEN 

C  -  FIND  RANGE  DF  GRID  PLANE  AND  ROUND  TO  WHOLE  C6L5II 

HMAX  =  HC  ♦  .5  P  HORZ 
HMIN  «  HC  -  .5  *  HORZ 
VMAX  «  VC  ♦  .5  P  VERT 
VMIN  •  VC  -  .5  P  VERT 

HMAX  *  SIGN( AINTC A3S(HMAX)  /  CEL5IZ  ♦  .5001)  *  CELSIZ,  HMAX) 

HMIN  =  SIGN(AINT(A3S(HMIN)  /  CELSIZ  ♦  .5001)  *  CELSIZ*  HMIN) 

VMAX  -  SICNCA1NTC ASS(VHAX)  /  CELSIZ  ♦  .5001)  *  CELSIZ,  VMAX) 

VMIN  »  SIGN(AINT(A9S(VMlN)  /  CELSIZ  ♦  .5001)  *  CELSIZ,  VMIN) 

ELSE 

WRITE  («,  6150) 

READ  (P,  5070,  £33*163 »  END-170)  HMIN,  HMAX,  VM1N,  VMAX 
ENOIF 

NriORZ  *  (HMAX  -  HMIN)  /  CELSIZ  ♦  1.00D1 
NVERT  -  ( VM  A  X  -  V  M I N )  /  CELSIZ  *  1.0001 
NCELL  «  NHDf. Z  P  WEST 

C  -  PRINT  PRIMARY  SLOCK  OF  HARDCOPY  OUTPUT 

WRITE  (P,  6160)  AZIM,  EL?V, 

1  TMlN(l)  *  MMCVT,  T M IN ( Z )  *  NmCVT,  TKlN(3)  *  MMCVT, 

2  TMAX(l)  P  MMCVT,  THAXC2)  P  MMCVT,  TMAXC3)  P  MMCVT, 

3  TCEN(l)  *  MMCVT,  TC?N(2)  *  MMCVT,  TCENC  3 )  *  MMCVT, 

A  TLEN(l)  *  MMCVT,  TL  EN ( 2 )  *  MMCVT,  TLENC  3)  *  MMCVT, 

5  SACK  *  MW'VT,  CELSIZ  P  MMCVT,  H3RZ  *  MMCVT, 

6  VERT  P  M “C  V  T  *  HC  «  MMCVT,  VC  P  MMCVT,  HMIN  P  MMCVT, 

7  HMAX  *  MMCVT,  V M I N  P  MMCVT,  VMAX  P  MMCVT,  NHO fi I , 

8  NVERT,  NCELL 
WRITE  (P,  a  1 7 0)  ISEED1 
IF  (RAYLOC  .EG.  C)  Then 

WRITE  (P,  6120 
ELSE  IF  CRAYLOC  .E3.  1)  THEN 
WRITE  (P,  6190  'RAYFILE  ' 

ELSE  IF  (RAYLOC  .ST.  1)  ThEN 
WRITE  (P,  6200) 

ENOIF 

IF  (PRTALL)  THEN 

WRITE  (P,  6210)  AZIM,  ELEV,  CELSIZ  P  MMCVT 
ENOIF 

WRITE  (1,  6210)  AZIM,  ELEV,  CELSIZ  P  MMCVT 

C  -  EITHER  READ  IN  RAY  PARAMETERS  DR  CALCULATE  THEM 

70  IF  (.NOT.  MX  R  AYS )  THEN 

READ  (51,  50  83,  ER3  =  180,  END-190)  HC  EN ,  VCEN,  H,  V,  IHIV,  EQV 
GCTC  100 
ENOIF 

IF  (ISPOT  .GT.  0)  THEN 
J  -  MOD( INPCEL,  4)  ♦  1 
IF  ( J  .EG.  1)  Then 

READ  (P,  5090,  ERR-160,  END-170)  (CH(Z),  CV(I),  I  -  1,  *) 
ENDIF 
H  «  CH( J) 

V  «  CV(J) 

HC  EN  *  SIGNCAINT(AES(H  /  CELSIZ)  ♦  .5)  *  CELSIZ,  H) 

VCEN  -  SIGNCAINT(ASS(V  /  CELSIZ)  *  .5)  P  CELSIZ,  V) 

WRITE  (P,  6220)  H  P  MHCV7,  V  *  MMCVT 
I H I  V  -  0 
ENOIF 
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non  nonnn 


VCEN  =  VM  AX 

C  -  PSEU00-LC3P  FOR  VERTIC4L  SCAN  -  SCANS  TCP  TO  30TTCH 

I  V:  RT  *  1 
80  HCEN  =  UMAX 

C  -  PSEJ03-L03P  FOR  HORIZONTAL  'CAN  -  SCANS  RIGHT  TO  LEFT 

IHGAZ  =•  1 

SO  1 r  CRAYL3C  .Si.  0)  THEN 

C -  CHOOSE  RANDOM  POINT  IN  CELL 

IRANV  *  10.  P  RANC'l) 

IRANI  “  10.  P  RANC-1) 

I H I V  s  10  P  IR4NH  ♦  IRANV 

V  =  VCEN  ♦  C  ELS  I Z  *  (.1  *  FLOAT (IRANV)  -  .45) 

M  *  hCEN  ♦  CELSIZ  p  C.l  *  FLOATCIRANH)  -  .AS) 

else 

c  -  choose  center  of  cell 

H  S  HCEN 

V  «  VCEN 
IMIV  =  0 

£  NO  I  F 

C  — *  JIGGL:  RAY 
100  CALL  TRQPICC-P) 

CAuL  CAOSSCdP,  HP,  W3> 

C  -  ROTATE  m,v  TO  COORDINATE  SYSTEM  Or  TARGET 

X£(  1 )  =  -V  u  CAZIM  *  SELEV  -  H  *  SAZIM  ♦  SJIGLE  *  rf?(I)  - 
1  SAC*  P  W3C1) 

Xa(2)  =  -V  P  SAZIM  a  SELEV  ♦M*  CAZIM  ♦  SJIGLE  P  WP(2)  - 
1  EACH  «  ¥3(2) 

X  3  (  3  )  “VP  CELEV  ♦  SJIGLE  P  «P(3)  -  BACK  P  N3(3) 

TRACK  CENTER  RAY 

DISTANCE  ALONj  SMCTLINE  TC  EIRST  CONTACT 

CRAYLEN  .GT.  0)  «  =  >  FIRST  CONTACT  IS  3ETUEEN  START  POINT 

A  NO  GRIG  PLANE 

CRAYLEN  .LT.  0)  «=)  FIRST  CONTACT  IS  BEYOND  GRID  PLANE 
CALL  FIRST 
RAYtENCOi  0)  “  0. 

ITHmITCC,  0)  “  -1. 

IF  (IREGON  .GT.  0)  THEN 

RAYLENCO,  0)  *  EACK  -  OFIRST 

CALL  UN2CLIRF0  •»  IREGON  -  1,  ITMHITCO.  0),  IOUM) 

£  N  o :  F 

TRACK  PLANETARY  RAYS 

NHhITS  IS  THE  NUM3ER  OF  PLANETARY  HAYS. THAT  HAVE  INTER  SECT  £0 
THE  TARGET 
NHHITS  s  0 

ANGINC  «  6.283135  /  FlC ATCNPLNTS ) 

00  120  PLANET  »  1,  NPLNTS 

SANGLE  -  (PLANET  -  1 >  P  ANGINC 
DC  110  CYL  *  1.  NMCYLS 

PlnETh(CTL,  PLANET)  «  H  ♦  CRAOCCYL)  P  SINCSANGLE) 
PLNETVfCYL,  PLANET)  ■  V  ■*  CRAOCCYL)  P  COSCSANGLE) 

XSC1)  *  -PLNETVCCYL.  PLANET)  *  CAIIH  «  SELEV  - 

1  PLNETHCCYL.  PLANET)  *  SAZIM  ♦  SJISLS  *  WPC1)  - 

2  BACK  P  W3( 1 ) 

X  3C  2 )  «  -PLNETVCCYL.  PLANET)  P  SAIIM  P  SELEV  ♦ 

1  PLNETHCCYL.  PLANET)  P  CAZIM  ♦  SJIGLE  P  WpC2)  - 

2  BACK  P  W3C2) 

X3C3)  «  PLNETVCCYL,  PLANET)  P  CELEV  ♦  SJIGLE  P  HP(3)  - 
1  BACK  P  m3 (3) 


31 


no 

120 


1 

2 


C 


1 

2 

3 


1 

2 

3 

130 

140 


c 


CALL  FIRST 

RATLENCCYL.  HAN'T)  .  0. 

ITMHITCCTL,  PLANET)  =  -1. 

IF  (IREGON  .GT.  0)  THEN 
NMHITS  -  NMHITS  ♦  1 

RAYlENCCYL,  PLANET)  -  SACK  -  DFIRST 

CALL  UN2CLIRF0  ♦  I R  EGON  -  1.  ITMHITCCTL,  PLANET),  IDUH) 
ENOIF 
CCNTINJE 
CONTINUE 

IF  CNHHITS  .Nr.  0  .OR.  ITMHITCO,  0)  .G2.  0)  THEN 
IF  CBRTALL)  THEN 

WRITE  CP,  t>230)  HCEN  *  HHCVT,  VC  E  N  *  HHCVT,  IHIV,  H  A  HHCVT, 
V  A  HMCVT,  RAYLENCO,  0)  *  HHCVT, 

ITMHITCO,  0) 

EN  OIF 


WRITE  (1,  6230) 
OUTPUT  HAIn-RAT 

a:  i a c  cyl  =  i, 

00  1 3 C  PLANET 
IF  C  0  R  T  ALL ) 
WRITE  CP, 


HCEN  *  HHCVT,  VCEN  *  HHCVT,  IHIV,  H  *  HHCVT, 

V  P  HHCVT,  RAYLENCO*  0)  *  HHCVT ,  ITMHITCO,  0) 

DATA 

NMCYLS 

=  1,  NPlNTS 

Then 

6240)  PlNEThCCYL,  PLANET)  *  MHCVT, 

PLNET VC  CYL,  PLANET)  *  HHCVT, 

RAYLENCCYL,  PLANET)  P  HHCVT, 

ITMHITCCTL,  PLANET) 


E  N  01  * 

WRITE  Cl.  6240)  PLNETHCCYL,  PLANET)  *  MMCvT, 

PLNETVCCYL,  PLANET)  P  HHCVT, 

RAYLENCCTL,  PLANET)  P  HHCVT, 

ITMHITCCTL,  PLANET) 


CONTINUE 


CQNTINJE 


ENOIF 

cNO  OF  CELL 


INPCEL  *  INPCEL  ♦  1 
IF  CISPCT  .GT.  INPCEL)  THEN 
GOTO  70 
ENOIF 

IF  CMARAYS)  then 

IF  CISPOT  .LE.  0)  THEN 

C -  NEXT  HORIZONTAL 

HCEN  »  HCEN  -  CELSIZ 
IHORZ  =  IHCRZ  •*  1 
IF  CIHCRZ  .LE.  NHORZ)  THEN 
GOTO  90 
ENOIF 

C  -  NEXT  VERTICAL  AT  ENO  CF  HORIZONTAL  SCAN 

VCEN  *  VCEN  -  CELSIZ 
IVERT  »  IVtfiT  4  1 
IF  CIVERT  .Lc.  NVERT)  THEN 
GOTO  BO 
ENOIF 
ENOIF 
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ELSE 

READ  (SI,  5030,  ERR*1S0,  EN0=19O)  HCEN,  VCEN,  H,  V,  IhIV,  EDV 
IF  (ECV  .NE.  999.9)  TMEN 
GOTO  100 
ENDIF 
ENDIF 

C -  END  OF  VIEW 

IF  (PRTALL  .AND.  NHASPS  .<JT.  1)  THEN 
WRITE  (P,  6250) 

EM  DIF 

WRITE  (P,  6260)  ASPECT 
WRITE  (F,  6  2  TO )  ISEEOl,  I  SEED 
ENDFILE  1 

IF  (IDLKP  .NE.  0)  THEN 
CALL  SECLH 
ENDIF 

C -  TIME  FOR  This  VIEW 

CALL  CLDCKS(ETIME) 

WRITE  (P,  62  30)  ASPECT,  ETIME  -  V T I H E 
150  CONTINUE 

C  -  END  OF  ALL  VIEWS 

END  FILE  1 
R  z  *  I N  D  1 

C  -  TCTAl  TIME  FOR  CYLINDER 

CALL  CLDCKSCzTIME) 

WRITE  (*,  6290)  ETIME  -  STIME 
RETURN 

C  -  HANDLE  READ  ERRORS 

160  WRITE  (P,  6300)  'INPUT' 

ST  D  P 

170  WRITE  (*,  6310)  'INPUT* 

STOP 

180  WRITE  (P,  63C3)  'RAYFIlE' 

STOP 

ISO  WRITE  CP,  6310)  'RAYfltc' 

STOP 

5010  FORMAT  (315) 

5020  FCRMAT  (A20,  215,  2=T.2,  2F9.3,  IS) 

5030  FORMAT  (3FI0.2) 

5  0*0  FDRMA7  (*=5.3) 

5050  FORMAT  (3=10.0,  315) 

5060  FORMAT  (3=5.0) 

5070  FORMAT  (4F10.0) 

5080  FORMAT  {/  10X,  *F10.*,  I*,  SX,  F6.1) 

5090  FORMAT  (8F10.0) 

6010  FORMAT  (’  NMASPS  *«UST  £  DUAL  1...'/ 

1  *  SU3R0UTINE  "SHOTCYL**  A90RTE0.') 

6020  FORMAT  ('0CPT10N  SET  TO  SUPPRESS  PRINTER  OUTPUT*) 

6030  FORMAT  ('  F21MXC  {',  13.  ')  >  TOTAL  NO.  OF  CYLS  (*,  13,  ')'/ 

i  *  subroutine  "smotcyl**  aborted.*) 

60*0  FORMAT  ('  ONE  OR  MORE  CYLINDER(S)  MUST  BE  SPECIFIED*/ 

1  *  IN  NON-DECREASING  ORDER...*/ 

2  *  SUBROUTINE  "SHOTCYL"  ABORTEO.*) 

6050  FORMAT  ('  *,  15,  A10,  A60) 

6060  FORMAT  C*  *,  L2,  A20,  2F7.2,  2F3.3,  15) 

60T0  FORMAT  (*  *,  *F8.3) 
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6060  FORMAT  <*0 -  GEOMETRY  OF  PROJECTILE  SKIN  (HH)  -- — ') 

6090  FORMAT  ('  RADIUS  CF  CYLINDER  ',  12,  '  '»  FIO. 2/ 

1  *  STEP3ACK  OF  CYLINDER  12,  *  FIO. 2) 

6100  FORMAT  ('  LENGTH  CF  NEEDLE-NOSE  PROBE  F10.2) 


6110 

FORMAT 

{'OPLANETARY  RAYS  PER 

CYLINDER  110) 

6120 

FORMAT 

('  ',  2110) 

6130 

format 

('  F10.2) 

61*0 

format 

('OTAROET  DESCRIPTION 

SAYS 

A 

=',  FIO. 3/ 

1 

V 

£•',  FIO. 3/ 

2 

m 

CELSIZ 

FIO. 3/ 

3 

'  BUT  RAYFIlE  SAYS 

k 

=',  FIO. 3/ 

4 

w 

£ 

FIO. 3/ 

5 

* 

CELSIZ 

=',  FIO. 3/ 

6 

*  subroutine  "smdtcyl" 

A£ 

ORTEO.') 

6190 

FORMAT 

C  *  0  0  BT ION  SET  TO  LIMIT 

GRID 

PLANE') 

6160 

FORMAT 

C'OAIIMJTH  ', 

FIO 

.3 

,  '  DEGREES*/ 

1 

'  ELEVATION  '♦ 

F  1  0 

•  3 

,  *  DEGREES*/ 

2 

'0 -  TARGET  - 

- 

X',  9  X ,  'Y*, 

3 

'  TARGET  MINIMUM  CMH) 

”  V 

3F10.3/ 

A 

'  TARGET  MAXI-JM  (MM) 

3^10.3/ 

9 

'  TARGET  CENTER  (MM) 

#  i 

3F10.3/ 

6 

'  TARGET  DIMENSIONS  (MM) 

• 

3F10.3/ 

7 

'0 - GRID  PLANE - 

-'/ 

6 

'  SACK  CFF  DISTANCE 

f  * 

FIO. 3,  '  MM'/ 

I 

'  CELL  SHE 

*  • 

FIO. 3,  '  MM'/ 

1 

'  HORIZONTAL  LENGTH 

F10.3,  '  MM'/ 

2 

'  VERTICAL  LENGTH 

'  « 

FIO. 3,  '  MM'/ 

3 

'  CENTER  (MM) 

#  ♦ 

2F10.3/ 

4 

'  HORIZONTAL  RANGE  (MM) 

• 

2F10.3/ 

5 

'  VERTICAL  RANGE  (MM) 

« 

2F10.3/ 

6 

'CNUM3ER  hCRZ  cells 

#  t 

110/ 

7 

'  NUMBER  VERT  CELLS 

9 

110/ 

S 

'  NUMBER  OF  CELLS 

9  % 

110) 

6170  FORMAT  ('OFIRST  SEED  CCS  RANDOM  NUMBER  GENE  RATER  * ,  I12> 
6160  FORMATCOOPTION  SET  TO  COMPUTE  RANDOM  POINT  IN  CELL') 

6190  FDRMATC 'OCPTION  SET  TO  READ  R^JNTS  FROM  FILE  A,  "") 

6200  FORMAT  C 'CCPTI ON  SET  TO  CHDDSE  CENTER  OF  CELL') 

6210  FORMAT  ('O',  3F10.2) 

6220  FORMAT  ('OS9ECIFIED  HCSI-*,  F10.2,  '  VERT-',  F10.2) 

6230  FORMAT  (’O',  2FB.2,  13,  3F8.2,  16) 

62*0  FORMAT  ('  *,  3FS.2 ,  16) 

6290  FORMAT  ('  999 . 9 ' , Ti X , ' E ND' > 

6260  FORMAT  C'lENO  OF  CASE',  IS) 

6270  FORMAT  C'OFIRST  SEED  FOR  RANDOM  NUMBER  GENERATOR',  112/ 

1  '  NEXT  SEEO  FOR  RANDOM  NUMBER  GENERATOR',  112) 

6260  FORMAT  C'OTIHE  FDR  CASE',  IS,  F9.3,  '  SECONDS'//) 

6290  FORMAT  C'OTOTAL  TIME  FOR  CYLINDER',  F9.3,  '  SECONDS') 

6300  FORMAT  ('OREAD  ERROR  ON  FILE  ',  A / 

1  *  SUBROUTINE  "SHOTCTL"  ABORTED. ') 

6310  FORMAT  ('OSOF  ENCGJNT EREO  ON  FILE  A/ 

1  '  SUBROUTINE  "SMDTCYL"  ABORTEO.') 

END 
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SUBROUTINE  SEEKVEW  CWH:NC£*  FNAME,  AZIMUTH) 


I  N  T  E  G  £  3 
CHARACTER  #7 
REAL 
LOGICAL 


WHENCE 
FnAmS 
AZIMUTH, 
EX  I  ST  F 


EOV,  VAZZM 


IN3UIHE  (F;l£=FVAME,  EXIST=EXISTF) 
IF  (.NOT .  EXISTF)  THEN 
WRITE  {*,  610)  FNAMj 
STOP 
ENG  I  F 

VAZIM  =  -1 
10  CLOSE  (WHENCE) 

OPEN  (WHENCE,  fil£*FnamE) 

READ  (kHENCE,  510,  END*30)  VAIIH 
IF  (VAZIM  .EO.  All  HUT  H )  THEN 
BACKSPACE  kH.EvCE 
RETURN 
ENCIF 

20  REAS  (WHENCE,  52C,  £ND*4C)  EOV 
if  (Eov  .EC.  959.S)  then 
R E a G  (whence,  *,  ESD=10) 

WRITE  (*,  620)  AZIMUTH,  FNAMr 
STOP 
E  NO  IF 
GOTO  20 

30  IF  (VAIIH  . SE .  0)  THEN 
WRITE  (*,  630)  VAZIM 
ELSE 

WRITE  («=,  60) 

£  SO  I F 

WRITE  (Ft  650)  AZIMUTH,  WHENCE 
STOP 

A 3  IF  (VAIIH  .iE.  0)  THEN 
WRITE  (*,  633)  VAZIM 
ELSE 

WRITE  (6,  6  A  3 ) 

E.N3IF 

WRITE  (*,  660)  AZIMUTH,  WHENCE 
STOP 


510  FORMAT 
520  FORMAT 
610  FORMAT 
1 

620  FORMAT 
1 
2 

630  FORMAT 
6* 0  FORMAT 
650  FORMAT 
1 

660  FORMAT 
1 
2 

END 


CF5-0) 

(/  6  2  X ,  F6.1) 

('OFILE  A,  NOT  FOUNO...'/ 

'  SUBROUTINE  "SEEKVEW"  A90RTED') 

C'ONECESSARY  *EOF  WAS  NOT  FOUND...'/ 

'  DIO  NOT  FIND  F5.1,  '  DEGREE  VIEW  ON  FILE 
'  SUBROUTINE  "SEEKVEW"  ABORTED') 

COLAST  VIEW  READ  WAS  F5.1,  '  DEGREES') 
('OREAD  NC  VIEWS') 

('  DID  NOT  FIND  ',  F5.1,  ’  DEGREE  VIEW  ON  UNIT' 
'  SUBROUTINE  "SEEKVEW"  ABORTED') 

('  NECESSARY  999.9  WAS  NDT  FOUND...'/ 

*  0 ID  NOT  FIND  ',  F5.I,  '  DEGREE  VIEW  DN  FILE 
'  SUBROUTINE  "SEEKVEW"  ABORTED') 


*  6/ 

15/ 

*  A/ 
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PROGRAM  NISCIR 
C 

C  COMPUTES  EFFFCTIVc  STANDOFF  CELL-BT-CELL  FOR  A  PARTICULAR 

C  THREAT/T ARSET  COMBINATION.  uses  output  fro*  shotctl  to 

C  DETERMINE  intersects  for  CYLINDER(S)  representing  threat. 

c 

C  WRITTEN  3T  PAUL  TANEN2AU**  IN  WHQH  ALL  BLAHE  LIES 

c  attn:  slcbr-vl-s 

C  APS,  HO  21005-5066 

C 

INTEGER  CELL,  CTL ,  CTLI,  FNCHAR,  FZ1MXC,  FI2HIT,  IKAIN,  IPLNET, 

1  NHASPS,  NMCELLS,  NHCYLS,  NMHITS,  NPLNTS,  PLANET,  READNO 

REAL  All*,  CO,  CELLPH,  CELSII,  CSTPBA*  DELAY,  DIST,  DLYOIS,  OHAIN, 

1  OPLNET,  ELEV,  GAP,  MCEN,  POIST,  PGAP,  SIGX,  SIGY,  STOOFF, 

2  VCEN,  VEL,  A,  XOC,  Y,  YOC 
LOGICAL  DEJUG,  EXISTF,  PPTALL 
CHARACTER  67  'NAME 
CHARACTER  610  OAT 

CHARACTER  620  RUNNAM 
CHARACTER  660  ITITLc 

Character  65a  shotln 

COMM!?*  /DEEUG/  OE3JS 
C  JMMC  VI  /ASPECT/  AZIN,  ELEV 
CGM-ICN  /CELSI :/  CfLSIZ 
CGKMCN  /0ISBCG/  XOC,  YOC 
COMMON  /OCSPER/  SISX,  SIGY 
COMMON  /VMS/  NMCELLS,  NHCYLS.  NPLNTS 
CCv«CVi  /RDUNC/  CD,  STDCFF 
OlHENSIJN  CST93K(3:iO) 

C  -  BEGIN  EXECUTION  - 

C -  C  a?N  FILES 

-  I?  THE  DATA  FILE  DDES  NOT  EXIST  THEN  ABORT 

(ANY  TRAILING  BLANKS  IN  REPRESENTATION  OF  FILE  NAME 

-  MUST  2E  IGNORED) 

FNAHE-  'TAPtl  ' 

CO  10  FNCHAR  «  1,  6 

IF  (FNAMECFNCHAR  ♦  1JFNCHAR  +  I)  .ED.  '  *)  THEN 
GOTO  20 
ENDI  F 

10  CONTINUE 

20  IN  GUI  RE  (FILE«FNAME(:FnCHAR),  E XI ST-EX I STF) 

IF  (.NOT.  EXISTF)  THEN 
WRITE  (6,  6310)  FNAME 
STOP 
ELSE 

OPEN  (SI,  FILE=FNA*E( JFNCHA*)) 

REWIND  51 
ENOIF 

CPE N  (61,  FlLE-'XXXSTRT') 

OPEN  (62,  FIL£=*XXXWTD*) 

DPEN  (63,  FILE='RESULT') 
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C  -  READ  HEADER  LINE  FOR  THIS  RUN 

READND  =  1 

READ  (51,  5C10,  ERR*110,  END=120)  NMASPS,  OAT,  XT  IT  LE 
WRITE  (63.  4020  NMASPS,  OAT,  ITITLE 
IF  (PRTALL)  THEN 

WRITE  (A,  6020  >  NMASPS,  DAT,  ITITLE 
ENDIF 

C  -  READ  PROJECTILE  DATA 

READND  ~  2 

READ  (  51,  5020,  ERR=90.  END*100)  PRTALL,  R'JNNAM,  CD,  STDCFF,  VEL, 
I  OELAT,  FZlMXC 

DEoJG  *  .FALSE. 

IF  ( ( V  £L  *  DELAT  A  FZlMXC  .EC.  0.0  .  AND. 

I  (VEL  ■*  LELAT  ♦  FZlMXC  .N-.  0.00  THEN 

WRITE  (A,  6030 
STOP 
ENDIF 

READND  =  3 

READ  (  31  ,  5330,  ERR  =  9C,  END=1CQ)  SIGX,  SIGY,  XDC,  T  DC 
IF  ( D E  5 L O  Then 

■  RITE  (A,  6040)  RJNNAM 
ENDIF 


C  -  READ  PROJECTILE  GEOMETRY  INFORMATION 

READND  =  4 

READ  (  51,  5040  ,  ERR >1 1 D,  END* 1 2  D )  NMCYLS,  NPLNTS 
READND  *  5 

DC  30  CTL  *  1,  NMCYLS 

READ  (51,  5050,  ERR*110,  END*120)  CSTPEX(CTL) 

30  CONTINUE 
READND  =  t 

READ  (51,  5050,  ERS*110,  END=120)  CSTPSK(O) 

C  -  READ  HEADER  LINE  FDR  THIS  VIEW 

READND  =  7 

READ  (  51,  5  06  C ,  ERR*110,  END  =  120)  aIIm,  ELEV,  CELSIZ 

WRITE  (63,  6050)  AZIM,  ELEV,  CELSIZ,  FZlMXC 
IF  (PRTALL)  THEN 

WRITE  (4,  6050)  AZIN,  ELEV,  CELSIZ,  FZ1NXC 
ENDIF 

C  -  FIND  THE  AINPOINT  (ASSJMSD  TO  5E  CENTER  DF  PRESENTED  AREA) 

CALL  CPA  (FNAHEOFNCHAR)) 

C  -  PRINT  OUT  RUN  INFORMATION 

WRITE  (A,  «063)  ITITLE,  CAT,  AZIM,  ELEV,  CELSIZ,  RUNNAN,  CD, 

1  STDOFr,  VEL,  OELAT,  SIGX,  SIGY,  ROC,  TOC,  NMCYLS, 

2  FI1HXC,  NPLNTS 
00  40  CTL  »  1,  NHCTLS 

WRITE  (A,  6070)  CTL,  CSTP3K(CYL> 

40  CONTINUE 

IF  (CSTP3KI0)  . GT.  0)  THEN 
WRITE  (a,  60  BO )  CSTP3X(0) 

ENDIF 

WRITE  (6,  5090) 
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oooooo  O  o  o 


C  -  DETERMINE  HDw  FAR  PROJECTILE  WILL  HDVE  DUPING  FUZE  DELAY 

C  -  (DISTANCE  ■  HATE  *  TIME) 

DLT01S  «  V-L  *  DELAT 

C  -  READ  IN  MAIN  SrtDTLlNE  eCR  NEXT  CELL 

CELL  -  0 
50  READNO  »  0 

READ  (51,  5070,  £RR=110,  END=SO)  SH2TLN 
READNO  *  9 

READ  (SHCTLN,  5030,  ERR-130,  END-140)  HCEN,  VCEN,  X,  Y,  CHAIN, 

1  IMAIN 

IF  (IMAIN  .EC.  -1)  THEN 
NHHITS  =  0 
ELSE 

NMhITS  »  I 
ENDIF 
GAP  =  0.0 
CYL1  «  0 
FZZhIT  *  0 

-  READ  IN  SHDTLINES  FOR  RAYS  CONSTITUTING  CYLINDER(S)  AND 

DETERMINE  WHICH  PLANET  RAY,  IF  ANY,  REPRESENTS  THE 
•—  REGION  OF  THE  PROJECTILE  WHICH  FIRST  IMPACTS  THE  TARGET 
DC  70  CtL  =  1,  NMCYLS 

DO  60  PLANET  =  1,  NPLNTS 
READNO  *  10 

READ  (51,  5  C  7  0 ,  ERR-110,  END  =  120)  SHOTLN 
READNO  *  11 

READ  CSHDTLN,  509C,  ESR-130,  END-140)  DPLNE7,  IPLNET 
IF  (IPLNET  .EC.  -1)  THEN 
GCTD  60 
ENDIF 

NHHITS  *  NHHITS  ♦  1 

PGA F  *  CPLNET  -  OMAIN  -  CSTP3XCCYL)  -  CSTPSK(O) 

-  (PGAP  .GT.  0)  *«>  THIS  PLANET  HITS  THE  TARGET  'BEFORE*  THE 

CENTER  OF  THE  PROJECTILE  DOES. 

(FZZHIT  .EC.  1)  =  ->  THE  PRIMARY  FUZE  HITS  THE  TARGET  FIRST 

CI.E.  0  <-  CYL1  <*  FUMXC),  BUT  THE 
SECONDARY  FUCE  HITS,  TOO,  BEFORE  THE  FUZE 

-  DELAY  IS  UP. 

IF  ((PG&P  .GT.  GAP)  .OR.  (NHHITS  .EC.  D)  THEN 
GAP  ■  PGAP 
CYLl  ■  CTL 

ELSE  IF  ((CYL  .GT.  FZ1MXC)  .AND.  (CYLl  .IE.  FZ1MXO)  THEN 
POIST  »  DHAIN  ♦  GAP  ♦  CSTPSK(O)  ♦  CSTPSKCCYL)  -  DPLNET 
IF  (POIST  .LT.  DLYDIS)  THEN 
FZ2HIT  «  1 
£  NO  IF 
ENDIF 

60  CONTINUE 

70  CONTINUE 

C  -  TURN  OFF  FZ2HIT  FLAG  IF  SECONDARY  FUZE  HITS  TARGET  FIRST 

IF  (CYLl  .GT.  F 2 1MXC )  THEN 
FZ2HIT  -  0 
ENDIF 
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C -  OEFlNE  ACTUAL  STANOCFF  I N  UMTS  OF  CHARGE  DIAMETERS 

CELL  =  CELL  ♦  1 
IF  (DMAIN  .EC.  0)  THEN 

C -  MARX  T  ri  2  S  CELL,  MEANING  THE  MAIN  SHDTlINE  MISSED  THE  TARGET 

DIST  =  0 
ELSE 

DIST  *  (GAP  ♦  ST  GOFF )  /  CD 
ENDIF 

C  -  DETERMINE  THE  PR03AaiLlTY  THAT  THIS  CELL  WILL  SE  MIT 

CELLPH  «  PHITCX,  Y) 

C  -  STORE  THIS  CELL'S  LOCATION  AND  CONTENTS  IN  FILE  'XXXSTRT', 

C  ITS  LOCATION  AND  WEIGHTED  CONTENTS  IN  FILE  'XXXWTO*,  AND 

C  -  ITS  CELL-LOCATION  AND  CONTENTS  IN  FILE  'RESULT' 

WRITE  (SI,  6100)  MCEN,  VCEN,  DIST,  1. 

WRITE  (62,  6100)  HCEN,  VCEN,  DIST,  CELLPM 

WRITE  (S3,  4110)  HCEN,  VCEN,  DIST,  CYL1,  FI2MIT,  CELLPH 

IF  (PRTALL)  THEN 

C  WRITE  (P,  6110)  HCEN,  VCEN,  DIST,  CYL1,  FI2HIT,  CELLPH 

WRITE  (P,  6110)  X,  Y,  DIST,  CYL1,  FZ2HIT,  CELLPH 
6NDIF 
GOTO  50 

c -  Close  files 

80  CLOSE  (51) 

CLOSE  (61) 

CLOSE  ( 62 ) 

WRITE  ( 63 ,  6120)  9)9.9 
CLOSE  (63) 

C  -  PRINT  HISTOGRAMS  OF  FIRINGS  VS.  ACTUAL  STANDOFF 

CALL  HlSTGG  ('XXXSTRT',  .FALSE.,  RUNNAM) 

CALL  n 1  ST  03  C'XXXwTD',  .TRUE.,  RUNNAM) 

C  -  PRINT  SILHOUETTE 

CALL  SILDET  ('XXXSTRT',  .TRUE.,  RUNNAM) 

IF  (DE6UG)  THEN 
WRITE  (P,  6130) 

ENDIF 

STOP 

C  -  HANDLE  READ  ERRORS 

90  WRITE  (*,  6140)  READND,  'ERROR  IN  FILE  'INPUT' 

STOP 

100  WRITE  (6,  6140)  REAONO,  'EOF  IN  FILS  ',  'INPUT' 

STOP 

110  WRITE  («,  6I4C)  READND,  'ERROR  IN  FILE  '*  FNAME 
STOP 

120  WRITE  («,  6140)  READND,  'EOF  IN  FILE  *,  FNAME 
STOP 

130  WRITE  (*,  6140)  READND,  'ERROR  IN  ',  'SHOTLN' 

STOP 

140  WRITE  (*,  6140)  REAONO,  'EOS  IN  ',  'ShDTLN' 

STCP 
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5010  FORMAT  (IX,  15,  A10,  A60) 

5020  FORMAT  (IX,  L2,  A20,  2F7.2,  2F6.3,  15) 

5030  FORMAT  (IX,  *F3.j) 

5  0*0  FORMAT  (IX,  2110) 

5050  FORMAT  (IX,  F10. 2) 

5060  FORMAT  (IX,  3F10.2) 

5070  FORMAT  (A) 

5080  FORMAT  (IX,  2F3.2,  3X,  3FS.2,  It) 

5090  FORMAT  (17X,  Ff.2,  16) 

6010  FORMAT  (  *1FIL  £  A,  *  NOT  FOUND*/*  PROGRAM  "MISFIR"  ABORTED*) 
6020  FORMAT  (IX,  15,  A10,  A60) 

6030  FORMAT  ('13AD  PRIMARY  FU2E  DATA*/'  PROGRAM  "MISFIR"  ABORTEO*) 
60*0  FORMAT  ('1ENTER  "MISFIR".  THE  TARGET  IS  *,  A) 

6050  FORMAT  (IX,  2F3.1,  2*X,  F8.1,  12) 


6060  FORMAT 

(*1 - PROGRAM  Ml  SfI 

R  — 

- / 

1 

'CTARGET  IS  *,  A60,  * 

RUN 

CN  *,  A10,  / 

2 

*  AZIMUTH 

*,  F6.2,  *  OEGRESJ 

3 

*  ELEVATION 

',  Fa. 2,  *  DEGREES 

* 

*  CELL  SHE 

',  F7.2,  *  MM*,  / 

5 

*o?r:ject:le  is  *,  A2o 

•  / 

6 

'  CHARGE  DIAMETER 

*,  F 6. 2  ,  '  MM',  / 

7 

'  3UILT-IN  STANDOFF 

',  F6.2,  *  MM',  / 

8 

*  IMPACT  VELOCITY 

F6.2,  *  M/S*,  / 

5 

*  FUZE  DELAY  TIM? 

',  F  6.2 ,  *  MS',  / 

e 

'  X  DIS»ECSIDN 

0 

•  F  6  •  3  • 

l 

*,  t  DISPERSION 

F8.3,  / 

2 

*  X  CORRECTION 

0 

,  F  9  •  3 , 

3 

*,  Y  CORRECTION 

F  8 • 3 ,  / 

* 

'  NUMBER  OF  CYLINDERS 

*,  12,  / 

5 

*  LAST  CYL  IN  PRIMARY 

fu:e 

12,  / 

6 

*  NUMBER  OF  RAYS  PER  CYL 

*,  12,  /) 

6070  FORMAT 

('  STEP3ACK  OF  CYLINDER 

« 

12,  2 X ,  F6.2) 

6030  FORMAT 

C  LENGTH  DF  O-WIOTM  PR 

DEE 

*,  F6.2) 

6050  FORMAT 

(* - 

- —  — 

- *) 

6100  FORMAT  (3FE-2,  FS.5) 

6110  FORMAT  (IX,  F6.0,  F7.1,  F3.*,  215,  F8.6) 
6120  FORMAT  (F6.1) 

6130  FORMAT  (*  NORMAL  EXIT  OF  "MISFIR*") 

61*0  FORMAT  ( *1 R  EA  0  **,  12,  *E  NC3UNT  ERED  *,  2A/ 
1  *  PROGAm  "MISFIR**  ABORTED*) 

END 
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SUBROUTINE  CPA  (FNAME) 


INTEGER  NMCELLS,  NMCYLS,  NPLNTS,  PLANET,  READNO 

REAL  DJMMT,  X,  X  A I M  ,  X9AR,  XOC.  T»  Y  A  I M  ,  YBAR,  YDC 

LOGICAL  DEEUG 

CHARACTER  *7  fNAME 

COMMCN  /AIM/  XA IM,  TAIN 

C CH MON  /DEBUG/  DE5UG 

COMMON  /DISPCO/  XOC,  YOC 

COMMON  /NHS/  NMCELLS,  NMCYLS,  NPLNTS 

C -  BEGIN  EXECUTION - - - - — ---- 

IF  CO  E  BUG)  THEN 
WRITE  C*,  610) 

ENOIF 
X  3AR  *  0 
YEAR  *  0 
NMCELLS  *  0 

C  -  FOR  EVERY  BUNGLE,  READ  COORJINATES  OF  CENTER 

13  REAuNO  =  1 

R  E  AO  (51,  510,  E  ND  =  3  0  )  X,  Y 

C  — -  ENSURE  THAT  REST  OF  DATA  FOR  CURRENT  BUNDLE  IS  CK 
DO  20  PLANET  *  1,  NPLNTS  *  NMCYLS 
REAONO  1  2 

READ  C51 ,  520,  ERR*50»  SND=60)  DUMMY 
20  CONTINUE 


C 


CCH»UTE  XEAS  AND  TSAR 


X  3  A  R  *  XsAR  ♦  X 
Y3AR  =  Y  3  A  R  ♦  Y 
NMCELLS  *  NMCELLS  ♦  1 
GOTO  10 

30  IF  CNMCELLS  .Ei.  C)  THEN 
WRITE  ($,  620)  'ERROR! 


NO  BUNDLE  DATA  LOADED.' 


WRITE  (*,  620)  'SUBROUTINE  "CPA"  ABORTED.' 


STOP 

ENOIF 

X3AR  3  XBAR  /  NMCELLS 
Y  B AS  «  YBAR  /  NMCELLS 


C  -  CALCULATE  THE  COORDINATES  OF  THE  AIMPQINT... 

C  -  AIHPOINT  IS  CENTER  Oe  PRESENTED  AREA  CORRECTED  FOR  DISPERSAL 

XAIH  B  XBAR  *  XOC 
Y AIM  «  YBAR  ♦  YOC 
IF  (DEBUG)  THEN 

WRITE  («,  630)  XBAR,  TSAR,  XAIM,  Y  AIM 
ENDIF 

C  -  REPOSITION  FILE  TO  ORIGINAL  POSITION 

REWIND  51 

DO  AO  I  «  1,  NMCYLS  *  6 

READ  (51,  *,  Efift*5Q,  ENO-60) 

60  CONTINUE 

RETURN 
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c 


-  HANDLE  file  ERRORS 

50  WRITE  <*,  620  )  'IRE  AD  READN3, 

1  '  ENCOUNTERED  AN  ERROR  IN  FILE  *,  FNAME, 

2  '  WHILE  ATTEMPTING  TO  REAO  CELL  NUMBER*. 


GOTO  70 

60  WRITE  (f,  620) 
1 
2 

70  WRITE  (*,  620 
STOP 


*1 R E AD  **,  REAONO. 

*  ENCOUNTERED  EOF  IN  FILE  PNAKE » 

*  WHILE  ATTEMPTING  TO  REAO  CELL  NUM3ER*. 

*  5  UsR 0 UT I NE  "CPA"  ABORTED.* 


NMCcLLS 


NMCELLS 


510  FORMAT  (IX,  19X,  2=3.2) 

520  FORMAT  (IX,  16X,  F9.2) 

610  FORMAT  (*  ENTER  SUBROUTINE  "CPA".') 

6t3  FORMAT  (A,  : 1 2 ,  2  A .  JA,  15) 

630  FORMAT  ('  X  3 A  R  =  *,  F10.A,  '  T  9 AR  «  *,  F10.A, 

1  *  X A I M  »  R10.A,  *  TAIM  *  F10.A) 

END 
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FUNCTION  PrllT  CX,  Y) 


8c 4 L  CELSIZ,  PhIT,  POFXY,  CCFZI,  QCFZ2 ,  GCFZ3,  3CFZ4,  SICX,  SIGY, 
1  X,  XAIM,  X  G  8 1  3 ,  T,  TAIN,  YGRI3 

LOGICAL  OioUG 
COMMON  /AIM/  XAIM,  Y  A  I  M 
COMMON  / j £ s  UG/  OE3UG 
COMMON  /3ISPER/  SIGX,  SIGY 
COMMON  /CELSIZ/  CELSIZ 


IF  CSIGX  «  SIGY  .EJ.  O)  THEN 

WAITE  C*.  o)  'ZERO  DISPERSION  NOT  ALLOWED,' 

write  c * ,  *>  'function  "pmit"  asorteo.' 

STOP 

ENDIF 

C  -  CALCULATE  THE  PS03A3LITY  OF  HITTING  TME  LOCATION  CX,  Y) 

XGRID  =  CElSIZ  /  2 

YG8ID  *  CELSIZ  /  2 

S  OF  1 1  =  0  =  NCC  X  -  XAIM  ♦  XGRID)  /  SIGX) 

C0FZ2  =  0=NC(X  -  XAI-  -  XGRID)  /  SIGX) 

3CFZ3  *  0  =  N  (  C  Y  -  Y  A I M  ♦  YG8I0)  /  SIGY) 

30F l A  =  3FSCCY  -  Y  A I M  -  YGRIO)  /  SIGY) 

PCFXY  *  CCOFZl  -  :3=Z0)  *  CC3eZ3  -  SSFZ4) 

PHIT  =  POFXY 
IF  C  D  E  3uG )  THEN 

WRITE  C9,  610)  X,  XAIM,  XGRIO,  SIGX 
WRITE  C*,  o CO)  Y,  Y A I M ,  YGRI3,  SIGY 
WRITE  C*.  630)  COFZI,  Q0FZ2 

WRITE  (9,  60)  CCFZJ,  COFZA 

WRITE  Cp,  650  )  P  3  =  X  Y 

ENOIF 
RETuRN 

6 1 C  FORMAT  C*  X*',  =10.4,  '  XAlM=',  FID. 4,  '  XGRID*',  F10.4, 

1  *  SIGX*',  F 1 0 • 4 ) 

620  FORMAT  C'  Y  =  ',  F10.4,  '  Y  A I M* ' ,  F10.4,  '  YGRIO*',  F13.4, 

1  '  SIGY*',  F10.4) 

630  FORMAT  C'  •  OF  Z 1 »  '  ,  =10.4,  '  C3CZ2»*.  FI0.4) 

640  FORMAT  C'  wGF  Z3*  *,  FI0.4,  '  C3=Z4=',  F10.4) 

650  FORMAT  C'  POFXY*',  F10.4) 

ENO 
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non 


c 


610 


REAL  FUNCTION  OFN  (  X  ) 

FROM  HASTINGS  APPROXIMATIONS  FDR  DIGITAL  COMPUTERS 
(SORROWED  FROM  WILSON'S  FILE  'CETANK' 

□STAINED:  9  Jun  83) 

REAL  A9S0FX,  F,  X 
LOGICAL  DEBUG 
COMMON  /DE3UG/  DE9UG 


F  =  0 

ABSOFX  =  A3S(X) 

IF  CA8S0FX  .LT.  5)  THEN 

F  «  (((((.  5333E-5  ft  ASSOFX  ■»  .493906E-4)  *  A3SDFX  ♦  .380C36E-A) 

1  6  ABSOFX  ♦  .3032776263)  *  ABSOFX  ♦  .  02  1 1  <-  1 00  61 )  *  ABSOFX 

2  ♦  .0*98673469)  *  63S3FX  ♦  1 
F  *  . 5  /  CFSS16) 

ENDIF 

IF  (X  . GE .  0)  THEN 
F  *  I  -  r 
ENDIF 
DrN  «  F 

IF  COESUG)  THEN 

WRITE  (*.  610)  X,  F 
ENDIF 
RETURN 


FORMAT  ('  OFN (  ' t  FIC.4,  ')  =  *,  FID.*) 
END 
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SUBRDUT INE  HISTOG  (FNAME,  WTO »  TARGNA) 

INTEGER  9  I  N  ,  CELL  ,  CCLO,  COLUMN,  I  LA  9  EL ,  LINE,  NPLNTS.  NMBINS, 
1  NMCELLS,  NMCYLS 

REAL  blNSIZ  ,  CD,  CFLTES,  FREQ  *  HMAX,  HSCALF,  MAXVAL,  RLABEL, 

1  STD3FF,  VALUE,  V»4X,  VSC4LF,  WT 

LOGICAL  DEBUG,  WTD 
CHARACTER  *7  FNAME 
CHARACTER  P20  TARGNA 
CHARACTER  «1I0  IMAGE 
COMMON  /3E5UG/  DEBUG 
COMMON  /NMS/  NMCELLS,  NMCTLS,  NPLNTS 
COMMON  /ROUND/  C  3 ,  ST  33  F  F 

DIMENSION  FReQC-UHO),  IMAGEC52),  VALUE  ( 1  OODO)  ,  WT(IOOOO) 
PARAMETER  (CCL0*11,  DFLT3S=0.1J 


c  - 

IF  (DEBUS)  TrtEN 
WRITE  (A,  6010) 

EN3IF 

c -  initialize  ofltbs,  fbejo,  imageo,  maxval,  nhbins,  and  vmax 

DO  10  BIN  =  -1»  100 
FR5KEIN)  =  0 
10  CONTINUE 

00  20  LINE  =  1,  52 
IMAGE  (LINE)  =  '  ' 

20  CONTINUE 
MAXVAL  =  C 
NMBINS  =  C 
VMAX  *  0 

C  -  READ  VALLEO  AND  «T()...  DETERMINE  THE  LARGEST  ELEMENT 

C -  ■  IN  VALUEC) 

OPEN  (5,  FlLE=eNAM£) 

R E n INO  5 

READ  (5,  5010,  E5?=150,  EN3*150) 

1  ( V  A LUE ( CELL ) ,  wT(CELL),  CELL-I,  NMCELLS) 

CLOSE  (5) 

0 D  30  CELL  =  1,  NMCELLS 

MAXViL  »  MAXCMAXVAL,  VALUE(CELL)) 

30  CONTINUE 

C  -  FOR  EACH  CELL,  ADO  WTCCELL)  TO  THE  BIN  ASSOCIATED  WITH 

C  -  VALUE (CELL) 

IF  (MAXVAL  .LS.  IOC  *  OFLTSS)  THEN 

C  -  THERE  ARE  NO  MORE  THAN  100  6IN5J  EACH  IS  OF  SIZE  OFLTBS 

BINSIZ  »  OFLTBS 
ELSE 

C  -  THERE  APE  100  BINS;  EACH  IS  OF  SIZE  MAXVAL  /  100 

BINSIZ  »  MAXVAL  /  ICO 
ENDIF 

DC  AO  CELL  *  1,  NMCELLS 

IF  (VALUE(CELL)  .ED.  0)  THEN 
BIN  ■  -1 
ELSE 

BIN  -  VALUECCELL)  /  BINSIZ 
ENDIF 
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FRE3C3IN)  ■  FRE3C3IN)  ♦  WUCELL) 

NM5INS  «  MAX  CNM3INS,  BIN) 

IF  COtaUS)  THEN 

WRIT'  (*•  6020)  CELL.  VALUECCELL),  WTCCELL).  BIN,  NM  3  INS 
E  NO  I F 

60  CONTINUE 

MM AX  e  NM3INS  *  3INS1Z 

C  -  SCALE  THE  CISTRI3UTICN  CURVE  TO  FIT  ON  THE  PAGE 

DO  50  3IN  =  -1,  N“3 1  NS 

V  HA  X  =  MAXCVMAX,  FREOCeiN)) 

50  CONTINUE 

IF  JHMAX  *  VM  AX  .53.  0)  THEN 
WRITE  C*.  6030) 

STOP 

SNOIF 

HSCALr  »  100  /  HMA  X 
VSCALF  s  55  /  VM  AX 

C  -  PLOT  HISTOGRAM  INTO  IMASEC) 

00  TO  SIN  =  0,  NHSINS 

COLUMN  *  COLO  ♦  1C1  -  COIN  ♦  .5)  *  BINSII)  e  HSCALF 
00  60  LINE  «  50,  51.5  -  VSCALF  *  FREGC&IN),  -  1 

Imageclinexccljmnjcolumn)  =  *e* 

60  CONTINUE 
70  CONTINUE 

COLUMN  *  COLO  ♦  ICS 

DO  80  LINE  «  50,  51.5  -  VSCALF  6  FRE3C-1),  -1 

imassclinc)cc:lumn:::lu“n)  =  •*' 

80  CONTINUE 


c 


-  ADO  AXES  TO  IMASEC) 

DO  90  LINE  s  I,  SI 

IF  CMCOCLINE,  10)  .EO.  1)  THEN 

IMASEUINEXCCLO  -  1JC0L0  -  1)  ■  •-* 
ELSE 

IMASECLINEKCCLO  -  1:CCL0  -  1)  «  *:* 
END  IF 

50  CONTINUE 

DO  100  COLUMN  =  CCLO  -  Z,  COLO  ♦  100 
IF  CM30CCCLUMN,  10)  ,E3.  1)  THEN 
IMAGE  C51)CCCLUMN:CDLUMN)  *  *1* 

ELSE 

image  C51)CCCLUMN:CCLUMN)  • 

ENDIF 

100  CONTINUE 

CDLUMN  *  COLO  ♦  105 
IMAGECSl)CCCLUMN:CCLUMN)  «  'E* 


C  -  LABEL  AXES  IN  IMASEC ) 

IF  CMAXVAL  •££•  50)  THEN 
DO  110  LINE  *  1,  61,  10 

ILA3EL  ■  C  51  -  LINE)  /  VSCALF 

WRITE  C  IMASECLINEKCCLO  -  7SCCL0  -  3),  6060)  ILABEL 
110  CONTINUE 
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ELSE 

OZ  120  LINE  »  1,  *1,  10 

RLAStL  =  (SI  -  LINE)  /  VSCALF 

WAITE  (IMADE(LlNE)CCGLO  -  lOSCCLO  -  3),  60  50  >  RLA6EL 
120  CONTINUE 
E  NO  I F . 

IHAi£C51)(CCLO-2!C:LO-n  *  'o:' 

IF  (BINSIZ  .EC.  DFLT5S)  THEN 

WRITE  CI*AGE(52)(CCL0-2:),  6060) 

ELSE 

WRITE  <IHAGE(52)CCCL0-*:>,  6070)  (UMAX  *  COLUMN  /  1O0. 

1  COLUMN  =  100.  20,  -20),  0 

t  NO  IF 

C  — -  WHEN  DEBUGGING,  PRINT  CONTENTS  OF  EACH  8IN 

IF  (DEBUS)  TnEN 
WRITE  (P,  SOSO ) 

00  130  BIN  ■  -1,  NMBInS 

wRITE  (p,  6050)  eiN,  FREDC3IN) 

130  CONTINUE 
£  NO  I F 

.  C -  PRINT  IMAGEC) 

WRITE  (*,  6100) 

IF  (.NOT.  WTO)  THEN 

WRITE  (f,  6110) 

ELSE 

WRITE  (6,  6120) 

ENOIF 

WRITE  (P,  6130) 

WRITE  (a,  51*0)  TARGNA 

WRITE  (P,  6150)  HH  A  X  /  100.,  HMAX  PCD  /  100. 

WRITE  (*,  6160)  VMix  /  50. 

DO ’ 1*0  LINE  a  1.  52 

■  RITE  (P,  e 1 7 0 )  IMAGE(LINE) 

1*0  CONTINUE 

WRITE  (p,  6030) 

RETURN 

C  -  HANDLE  READ  ERRORS 

150  WRITE  (P,  6130)  'READ  SCREW  U».  SUSROUTINE  "MISTOG"  ABORTED.' 
STOP 

5010  FORMAT  ( 1 6 X ,  *3.2,  FS.6) 
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6010  FORMAT  ('  ENTER  SUBROUTINE  * *HJ ST 3 & ' *. ' > 

6020  FORMAT  {'  VALUE(',  Ia,  ')  -  ',  F10.3,  '  WT  •  *,  F10.3, 

1  '  SIN  IS  13,  '  OUT  OF  ',  13) 

6030  FORMAT  ('HERO  Or  NOMINATOR  '/'PROCEDURE  "hISTOG"  ABORTED') 
60a0  FORMAT  CIS) 

6050  FORMAT  OPE:. 2) 

6060  FORMAT  C'10',  9X,  '9',  9k,  '6',  9k,  '7',  9k,  '6'.  9k,  '3',  9k, 
1  'A',  9k,  '3',  9k,  '2',  9k,  '1',  9k,  '0') 

6070  FORMAT  (5CIPE3.Z,  12X),  AX,  II,  AX,  'M') 

6080  FORMAT  C'l') 

8090  FORMAT  ('  'FRJC(',  13,  ')  *  *,  FT. 2) 

6100  fORMAT  C'l',  SAX,  'HISTOGRAM') 

6110  FORMAT  (*  Six,  'NJM3ER  OF  CELLS') 

6120  FORMAT  C'  A2X,  'NJ“sER  OF  CELLS  (WEIGHTED  FOR  PH)') 

6120  FORMAT  (*  ',  ATX,  'VS  ACTUAL  STANDOFF  (CD)') 

6 1 A  0  FORMAT  {'O',  A  2  X ,  'TARGET  IS  '",  A,  "") 

6150  FORMAT  ('  ',  30X,  'HORIIONTAl  UMT  LENGTH  IS  ',  1PE9.2, 

1  'CD  (  *  ',  1PE9.2,  '  MM)') 

6160  FORMAT  ('  *,  a  2  X ,  'VERTICAL  UNIT  LENGTH  IS  ',  1PE9.2) 

6170  FORMAT  (11X,  A) 

fcleO  FORMAT  (*  ',  A) 

END 
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SUBROUTINE  SIL0E7  (FNAME,  M£7RIC»  TARJNA ) 


C  DERIVED  F  R  CM  WILSON'S  'SILPK*  SUBROUTINE, 

C  obtained.  a  ju ns  1933 

C  DRASTICALLY  REWORKED  BY  PAUL  T ANEN  BA’JM 

INTEGER  COL.  IH,  IV,  JVC,  LABEL,  RON,  STARTC,  StARTL.  STOPC, 

1  ST D PL  »  VALUE 

REAL  AZIM,  CD,  CEL  SI  I ,  ELEV,  PK ,  STDOFF ,  VAL,  NT,  X,  XMAX,  XMAX1. 
1  XMIN,  XW INI,  Y,  YHAX,  YMAX1,  YMIN,  YN  INI 

LOGICAL  DEBUG,  L  A  SpAG  ,  METRIC 

Character  r7  fname 

CHARACTER  <rZ  3LNK,  HLA5EL,  VLABEL,  HZERC,  IMAGE,  VZERO 
character  *20  targna 

DIMENSION  IMAGEC2S0,  130),  HL  ABEL  (250),  VLABEL(130),  VALO) 

LEVEL  2,  /BUFFER/ 

CCMmOn  /ASPECT/  AIIM,  ELEV 
COMMON  /SUFFER/  IMAGE 
COMMON  /OELSIZ/  celsii 
COMMON  /OESUG/  DEBUG 
COMMON  /ROUND/  CD,  STDOFF 
DATA  VZERO  /'--'/ 

DATA  HZERC  / . / 

DATA  E N  X  /  '  '/ 


IF  { 3£  3  US )  T  H  e  N 
WRIT;  (*,  610) 

SNDIF 

c -  3LANK-FILL  r!L  AS;  L  (  )  «  VLABELO,  ANC  IMAGEC) 

0010  COl  *  1,  250 
HLA3  EL  C  COL )  =  SLNK 
10  CONTINUE 

30  20  ROa  «  1,  133 

VLASEL(ROW)  *  BLNK 

20  CONTINUE 

DO  <.0  COL  -  1,  250 
DO  30  R 0 a  =  I,  130 

IMAGECCOL,  R  0  a )  *  3LNK 
30  CONTINUE 
«0  CONTINUE 

C  — -  DEFINE  DIMENSIONS  OF  OUTPUT 

XMAX  =  -125  *  CELSIZ 
XMIN  ■  125  *  CELSIZ 
YMAX  »  -56  *  CELSIZ 
YMIN  «  55  »  CELSI2 
XMAXI  •  XMIN 
XHINI  «  XMAX 
YHAXI  =  YMIN 
YMIM  a  YMAX 
LASPAG  ■  .FALSE. 
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C  -  READ  DATA  FOR  EACH  CELL 

OPEN  (5,  FILE'FNAHc) 

REWIND  5 

50  READ  (5.  51C,  £RR=90,  *ND*60)  X,  Y,  PR,  WT 
PR  ■  PR  *  WT 

C  -  REDEFINE  DIMENSIONS  OF  OUTPUT,  IF  NECESSARY 

XHAX  *  HAXCXMAX,  X) 

XMIN  =  MINCXMIN,  X) 

YMAX  ■  MAXCYMAX,  Y) 

YMIN  ■  MINCYMIN,  Y) 

C -  WRITE  INTS  HLA3ELO  AND  VLABELO 

COL  »  ((X  -  XMINI)  /  CELSIZ)  1 
LABEL  =  ABS(X)  /  CELSIZ 

CALL  NJ.NSITCHL  ABrLCCDL),  LABEL,  METRIC) 

ROW  *  ((TmaxI  -  Y)  /  CELSIZ)  ♦  1 
LABEL  *  A9SC Y)  /  CELSIZ 

CALL  NJMPITCVLAbELCRCW),  LABEL,  METRIC) 

C  — -  WRITE  INTO  IMA  SEC ) 

VALUE  =  (PR  *  10)  ♦  .5 

CALL  NUMRITdMA&ECCSL,  ROW),  VALUE,  .TRUE.) 
GOTO  50 

C -  PRIM  HEADERS 

60  CLOSE  (5) 

WRITE  (#■  &ZQ)  T  A  K  SNA 

WRITE  (?,  630)  SLE  V ,  AZIM,  CELSIZ 

C  -  CENTER  OF  TARSET 

CCl  ■  (CO  -  XMIM)  /  CELSIZ)  ■*  1 
ROW  =  (CYKAXI  -  0)  /  CELSIZ)  ♦  1 
HLA3EL(CCL)  =  MZERO 
VLASEL(ROk)  *  VZERO 

C  -  SET  MDRZ  AND  VERT  SPACIN5  On  PACE 

IH  =  (XHAX  -  XMIN)  /  CELSIZ  ♦  1 
IV  *  (YMAX  -  YMIN)  /  CELSIZ  ♦  1 
IVC  *  IV 

IF  (IH  .LE.  63)  THIN 

IH  *  (63  -  Irt)  /  2  ♦  1 
L ASP AC  =  .TRUE. 

ELS* 

I H  =  1 
£  NO  IF 

IF  (IVC  .LE.  55)  THEN 
IV  *  (56  -  IV)  /  2 
ELSE 
IV  *  1 
ENDIF 

C  -  DETERMINE  THE  WINDOW  OF  IMAOEO  TO  PRINT 

STARTC  *  (XMIN  -  XMINI)  /  CELSIZ  +  1  -  IH 
STDPC  •  STARTC  ♦  63 

STARTL  »  CYMAXI  -  YMAX)  /  CELSI Z  ♦  1  -  IV 
IF  (IVC  .LT.  56)  THEN 
STOPL  “  STARTL  ♦  55 
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ELSE 

STOPL  s  (YMAXI  -  TWIN)  /  CELSII  ♦  1  -  IV 
STARTL  *  Mt  XC  ST  A  RTL  ,  1) 

ST  A  ST  C  ■  HAXCSTASTC,  1) 

ENDIP 

C -  PRINT  SILHOUSTT  E 

DO  70  ROW  «  STARTL,  STOPL 

WRITE  C».  6400  VLA6ELCR2W),  CIMACECCCL,  ROW). 

I  COL  *  STARTC,  STOPC) 

70  CONTINUE 

WRITE  CP i  640 )  VLA3ELC1),  C HL A  9  EL C COL )  ,  COL  ■  STARTC.  STOPC) 
WRITE  CP,  650) 

C  -  OIO  THE  SILmOUETTE  FIT  ON  ONE  PACE? 

IF  (LASPiC)  THEN 
STOP 
£  NO  I F 

C  -  PRINT  REMAINDER  OF  SILHOUETTE 

WRITE  CP.  660) 

WRITE  <P,  670) 

STARTC  =  STOPC  ♦  1 
STOPC  =  STARTC  *  60 
00  90  ROW  *  STARTL,  STO^l 

WRITE  Cp.  640)  CIMACECCCL,  ROW),  COL  =  STARTC,  STCPC), 

1  VLA2ELC50W) 

SC  CONTINUE 

WRITE  C*.  6*0)  ChlA3;LC::L).  COL  =  STARTC,  STOPC),  VLABELC1) 
RETURN 


C -  (i ANIL E  READ  ERRORS 

90  WRITE  C  ? »  650  )  *IR=A0  5CRE-UP...  SJ5R0UTINE  "SILOET"  AeORTEO.' 
STOP 

510  FORMAT  (3*5.2,  F5.6) 

610  FORMAT  C'  ENTER  SUaROjTINE  "SILOET**.') 

620  FCRMATC'l *,  4  8  X ,  'CELL-BY-CELL  STANDOFF  IN  TENTHS  OF  A  CO*,  /, 

1  S3X,  'TARCET  IS  ',  A20 ) 

630  FORMATC*  *,  29X,  Fj.O,  '  DEGREES  ELEVATION  F5.0, 

1  *  DECREES  AZIMUTH  CELL  SHE  «  FT.O,  '  MM*) 

640  FORMAT  C1X.65A2) 

650  FCRMATC130X, '. ') 

660  FORMATC'l') 

6T0  FDRMATC1X,/) 

690  FORMAT  (A) 

END 
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SUBROUTINE  numrit  (whither,  WHAT,  HOW) 

INTEGER  WHAT 
LOGICAL  HOW 
CHARACTER  *2  WHITHER 


C 


IF  (HOW)  THEN 

If  (whAT  ,G5.  1C1)  THEN 
WRITE  (WHITHER,  610) 

ELSE 

WRITE  (WHITHER,  620)  WHAT 
ENOIF 
ELSE 

IF  CwHAT  .GE.  100)  Th=N 
WRITE  (WHITHER,  620)  99 
ELSE 

WRITE  (WHETHER,  620)  MOOCWMAT,  25)  *  A 
endif 
ENOIF 
RETURN 

610  FORMAT  ('A:*:') 

620  FORMAT  (12.2) 

E  NO 
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P.RCORAN  FUZES 


C 

C  COMPUTES  PROiAJILITY  THAT  A  2-PUZED  PROJECTILE'S  SECONOART  FUZ E 

C  WILL  STRIKE  TARGETS 

C  A)  BEFORE  PRIMARY  FUZE 

C  B)  OJRINS  PRIMARY'S  DFLAY 

r  r  t 

C  WRITTEN  3  Y  PAUL  YtNENSAUM,  IN  WHOM  ALL  BLAME  LIES 

C  ATTN!  SLC3R-VL-5 

C  APS,  MO  21005-5066 


INTEGER  ASPECT,  CYL1,  OLY,  0LY9A0,  FUZ  E2,  FZ1MXC,  NHASPS, 

1  N^OELLS*  R  E  A  ON  0 

SEAL  All-,  CELLPH,  CELSI2,  CIST,  DLYP,  ELEV,  FUZEZP,  HCEN,  VCEN 
CHARACTER  *10  DAY 
CHARACTER  *4  0  I T I T  L  E 

C -  5ECIU  EXECUTION  - 

OPEN  C51,  FUE»'FUZINP'3 


REAO.NO  =  1 

READ  ( 51  ,  510,  E  R  S  =  30 ,  END  =  4.9)  NHASPS,  DAY,  ITITLE 

C  ----  PROCESS  NEXT  VIEW 

CO  ZD  ASPECT  =  1,  NHASPS 

C  INITIALIZE  THINGS  FOR  THIS  VIEW 

NHCELLS  =  0 
FUZE2  =  0 
FJZEEP  =  0 
CL  T  =  C 
DLYP  =  0 
REACnC  =  2 

READ  <  51  ,  520,  E«R=30,  E  NO*  A  0  0  AIIH,  ELSV,  CELSIZ,  FZ1HXC 

FOR  EACH  CELL...  READ  DATA  AND  COMPUTE  PS 03 ABILITY  THAT 
SECONDARY  FUZE  WILL  STRIKE  FIRST  AND  PR39ABILITT  THAT  IT 
WILL  HIT  CURING  PRIMARY-FUZE  DELAY. 

RcACNJ  «  3 

READ  C  5 1  ,  530  ,  ERR=30,  END=A0)  HCEN,  VCEN,  DIST,  CYL1,  OLYBAD, 

CELLPH 

IF  CHCEN  .NE.  959. 1)  THEN 
IF  CCTL1  .CT.  CZ 1  NX  C  3  THEN 
Fu  Z  £  2  *  FUZ  E  2  ♦  1 
FUZE2P  *  Fu:*2P  ♦  CELLPH 
£  NOI c 

IF  (DLY3AD  • E2-  13  THEN 
DlT  »  iLI  ♦  1 
DLYP  ■  OLTP  ♦  CELLPH 
ENDIF 

NMCELLS  «  NHCELLS  ♦  1 
GOTO  10 
ENDIF 


C - 

C 

C - 

1  0 

1 
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C  -  OUTPUT  RESULTS 

WRITE  (6,  610)  ITITLE,  DAT »  AZIH,  ELEV,  CELSIZ,  NMCELLS, 

1  FUZ  E2»  ReAL(FUZE2)  /  RE AL(NMCELLS) *  eUZE2P, 

2  DLT,  REALCOLT)  /  REAL(NHCELLS),  OLTP 
20  CONTINUE 

STCP 

C  -  HANDLE  INPUT  ERRORS 

30  IP  C RE  AON 3  .£  0.  1)  THEN 
WRITE  (»i  623)  REAONS 
ELSE  IF  (RSAONO  .E2.  2)  THEN 
WRITE  («i  633)  REAONOt  ASPECT 

else 

WRITE  CP,  6*0)  READN3,  ASPECT,  NMCELLS  ♦  1 
ENDIP 
STOP 

*0  IF  (REA3N0  .E5.  1)  THEN 
WRITE  <*»  650)  READN3 
ELSE  IF  (REA0N3  .EO.  2)  THEN 
WRITE  (*,  660)  REA3N3,  ASPECT 
ELSE 

WRITE  (*,  6T0 )  REACN3,  ASPECT,  NMCELLS  ♦  1 
E  NO  IF 
STOP 


510  FORMAT  (IX,  15,  AID,  A60) 

520  FORMAT  (IX,  2F3.1,  2*X,  F3.1,  12) 

530  FORMAT  (IX,  F6.0,  FT.l,  F9.4,  215,  FB.6) 

610  FORMAT  Cl - PRO  SRAM  FUZES - ',  / 

1  'CTARGST  IS  '",  A  6  0 ,  '"  RUN  ON  ',  A10,  / 

2  *  AZIMUTH  *,  F6.2,  '  DEGREES', 

3  *  ELEVATION  *,  F6.2,  *  DEGREES', 

*  '  CELL  SIZE  *,  FT. 2,  '  MM',  / 

5  '  NONEMPTY  CELLS  THIS  VIEW  15,  / 

6  'OCELLS  HIT  BY  SECONDARY  FUZE  ',  15,  / 

T  '  PROBABILITY  GIVEN  A  HIT  ',  FA. 3,  / 

3  *  PROBABILITY  GIVEN  A  SHOT  ',  FA. 3,  / 

9  'OCELLS  HIT  DURING  FUZE  DELAT  15,  / 

6  *  PROBABILITY  GIVEN  A  HIT  ',  FA. 3,  / 

1  *  PR03A3ZLITT  GIVEN  A  SHOT  ',  FA. 3,  / 

2  *  - ') 


620  FORMAT  CllNPUT  ERROR  ON  READ  NUMBER  II,  / 

1  'PROGRAM  "FUZES"  ABORTED') 

630  FORMAT  CilNPUT  ERROR  ON  READ  NUMBER  II,  '  ASPECT*',  13,  / 

1  'PROGRAM  "FUZES"  ABORTED') 

6 A  0  FORMAT  CllNPUT  ERROR  ON  RE  AO  NUMBER  ',  II,  '  ASPECT-*,  Z3, 

1  '  CELL-',  IA,  /  'PROGRAM  "FUZES"  ABORTED') 

6S0  FORMAT  ('IP  RE  MATURE  EOF  ON  READ  NUMBER  ',  II,  / 

1  'PROGRAM  "FUZES"  ABORTED') 

660  FORMAT  ClPR-MATURE  EOF  CN  REAS  NUMBER  ',  II,  '  ASPECT*',  13,  / 
1  'PROGRAM  "FUZES"  ABORTED*) 

670  FORMAT  CIPREMATUFE  EOF  ON  READ  NUMBER  *,  II,  '  ASPECT-',  13, 

1  '  CELL-',  IA,  /  'PROGRAM  "FUZES"  ABORTED') 

END 


/ 

/ 
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APPENDIX  B 
Sample  Output 
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This  page  Left  Intentionally  Blank 


ENTER  USE-?  ROUTINE  SRANDX 


C  P  T  I  D  N  SET  T  D  SUPPRESS  PRINTER  OUTPUT 


- GEZM 

RADIUS  OF 
STEP3ACK 
RADIUS  CF 
STEP5ACK 
RADIUS  Dr 
STE03ACX 
RADIUS  DF 
ST  E°5ACK 
RADIUS  OF 
STEPbACA 


ETRY  DF  ORDJECTIUi 
CYLINDER  1 
OF  CYLINDER  1 
CYLINDER  2 
Or  CYLINDER  2 
CYLINDER  3 
DF  CYLINDER  3 
CYLINDER  4 
DF  CYLINDER  4 
CYLINDER  5 
OF  CYLINDER  S 


S  XI N  (MM) 
14.14 
5.96 
20.00 
20.00 
21.00 
153.00 
57.45 
176.79 
75.00 
225.00 


planetary  rays  ®er  cylinder 


s 


AZIMUTH  .000  DEGREES 

ELEVATION  .000  DEGREES 


-  TARGET  - 

X 

Y 

Z 

TARGET  MINIMUM  (MM) 

-3673.000 

-1642.000 

.000 

TARGET  MAXIMUM  (MM) 

6200.000 

1642.000 

2831.000 

TARGET  CENTER  (MM) 

1243.500 

.300 

1415.509 

TARGET  DIMENSIONS  (MM) 

9373.000 

3294.000 

2331.003 

-  GRID  PLANE  - 

SACK  2 =F  DISTANCE 

6594. COO 

MM 

CELL  SIZE 

100.000 

MM 

mcrizdntal  length 

3234.000 

MM 

vertical  length 

2331.000 

MM 

CENTER  (MM) 

.  009 

1415.500 

MCRIZDNTAL  RANGE  (MM) 

-1600.000 

1600.000 

VERTICAL  RANGE  (MM) 

.000 

2600.000 

NJM5ER  MDRZ  CELLS 

33 

NUMBER  VERT  CELLS 

29 

NuMsER  DF  CELLS 

957 

FIRST  SEED  FDR  RANDOM  NUM3ER  GENERATOR  0 

OPTION  SET  TD  COMPUTE  RANDOM  POINT  IN  CELL 
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END  3  F  CASE  1 

FIRST  SEED  FOR  RANDOM  NUMbER  GENERATOR  0 

NEXT  SEED  FDR  RANDOM  NUM3ER  GENERATOR  21337SS5 

TIME  FOR  CASE  1  73.594  SECONDS 

TDTAL  TIME  FDR  CYLINDER  73.SC6  SECONDS 
LEAVE  USER  ROUTINE  3RANCX 
END  DF  RUN 

./ 
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PROGRAM  MI3FIR 


TARGET  IS 
AZIMUTH 

elevation 

CELL  SIZE 


T-S2A  TANK  DESCRIPTION  CGIFT5)  RUN  ON  09/23/35 

.00  DEGREES 
.00  DEGREES 
100.00  MH 


PSG JECTILE  IS  DEMONSTRATION  RUN 


CHARGE  DIAMETER 
aulLT-lN  STANDOFF 
IMPACT  VELOCITY 
FUZE  DELAY  TIME 
X  DISPERSION 
X  CORRECTION 
VJMSER  OF  CYLINDERS 
LAST  CYL  IN  PRIMARY  FUZ 
NJMEER  CF  RAYS  PER  CYL 


100.00  MM 
250.00  MM 
175.00  M/S 
.50  MS 

400.000,  Y  DISPERSION 

.000,  Y  CORRECTION 

5 
2 
3 


400.000 

.000 


STEP3ACX  OF  CYLINDER  1 
3TEP3ACK  OF  CYLINDER  2 
STEP3ACK  OF  CYLINDER  3 
3TEP5ACX  OF  CYLINDER  4 
3TEP3ACX  CF  CYLINDER  5 


5.  36 
20.00 
I  53. CO 
176.79 
225.00 
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TiiiSI  i;  *::N3SiT«!-Tns  *'J>* 

KtiL  UhU  wiHiTH  IS  S.<3  =  -i»l  CO  (  *  5.jJ«*0l  **> 

Vi«T|C»L  lih  I T  lihiTd  IS  *.0S‘*OO 


r 


•  « 


1/ 

u 

r 


» 
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»3£»0I  i.»jc*01  3.S0i*0I  2.33t*01  1.17C*0l 


fi  ISTCGR14 

zr  CeLLS  (WEIGHTED  FO®  *hj 
VS  ACTUAL  STA*«23cF  (CO) 
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93E  •  01  t.«4E»01  J.jOE'QI  E.JJEOl  l.J7c«01 


ctiL-iY-CEii  standoff  in  tenths  of  a  cd 
TAAJSt  IS  OtNJNSTRAlIaN  RUN 

DtiHE'S  ELEVATION  0.  Ot&ilEcS  -CINjTrt  CELL  SHE  »  100. 


* 


NNMNNNNN^NpANAMMrt^fN^MOOOOeueQQ  | 
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161)10  1  31*1 119070401  060  SOLO  30*01  '*01020  JOWS063  706091  a  1 112131*1  Sid 


PROGRAM  FUZES 


TARGET  IS  'T-6ZA  TANK  DESCRIPTION  (GIFTS)  RUN  ON 


AZIMUTH 

.00 

DEGREES 

ELEVATION 

.00 

DEGREES 

CELL  SIZE 

100. CO 

MM 

NONEMPTY  CELLS  THIS 

VI 

EW 

677 

CELLS  HIT  3 Y 

SECONDARY 

FUZE 

146 

probability 

GIVEN 

A 

HIT 

.216 

PROBABILITY 

GIVEN 

A 

SHOT 

.121 

CELLS  HIT  DLR 

ING  FUZ 

w 

DELAY 

55 

PROBABILITY 

GIVEN 

A 

HIT 

.031 

PROBABILITY 

GIVEN 

A 

SHOT 

.053 

01/23/35 
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DISTRIBUTION  LIST 


No.  of 

Copies  Organization 

12  Administrator 

Defense  Technical  Info  Center 
ATTN:  DTIC-DDA 
Cameron  Station 
Alexandria,  VA  22314-6145 

1  HQDA 

ATTN:  DAMA-WSW 
Washington,  DC  20310 

1  HQDA 

ATTN:  DAMA-ART-M 
Washington,  DC  20310 

1  Commander 

US  Army  Materiel  Command 
ATTN:  AMCDRA-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333-0001 

1  Commander 

US  Army  Materiel  Command 
ATTN:  AMCLDC 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

1  Commander 

US  Army  Materiel  Command 
ATTN:  AMCDE-DW 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 


No.  of 

Conies  Organization 

10  Central  Intelligence  Agency 
Office  of  Central  Reference 
Dissemination  Branch 
Rm.  GE-47,  HQS 
Washington,  DC  20502 

1  Central  Intelligence  Agency 
OSWR/OSD/GPWB 
Washington,  DC  20505 

1  Commander 
Armament  Research  & 

Development  Center 
USAAMCCOM 
ATTN:  SMCAR-TDC 
Dover,  NJ  07801 

2  Commander 
Armament  Research  & 

Development  Center 
USAAMCCOM 
ATTN:  SMCAR-TSS 
Dover,  NJ  07801 

1  Commander 

Armament  Research  & 

Development  Center 
USAAMCCOM 
ATTN:  SMCAR-LC 
Dover,  NJ  07801 
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1  Commander 

US  Army  Armament  Munitions 
&  Chemical  Command 
ATTN.  SMCAR-LEP-L 
Rock  Island,  IL  61200 

1  Commander 

US  Army  Armament  Munitions 
&  Chemical  Command 
ATTN:  SMCAR-ESP-L 
Rock  Island,  IL  61200 

1  Director 

Benet  Weapons  Laboratory 
Armament  Research  & 
Development  Center 
USAAMCCOM 
ATTN:  SMCAR-LCB-TL 
Watervliet,  NY  12180 

1  Commander 

US  Army  Aviation  Research 
&  Development  Command 
ATTN:  AMSAV-E 
4300  Goodfellow  Blvd 
St.  Louis,  MO  63120 

1  Director 

US  Army  Air  Mobility  Research 
&  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  04035 


No.  of 

Conies  Organization 

1  Commander 

US  Army  Communications- 
Electronics  Command 
ATTN:  AMSEL-ED 
Fort  Monmouth,  NJ  07703 

1  Commander 

ERADCOM  Technical  Library 
ATTN:  DELSD-L  (Reports  Section) 
Fort  Monmouth,  NJ  07703-5301 

1  Commander 

US  Army  Missile  Command 
Research,  Development,  & 
Engineering  Center 
ATTN:  AMSMI-RD 
Redstone  Arsenal,  AL  35808 

1  Director 

US  Army  Missile  &  Space 
Intelligence  Center 
ATTN:  AIAMS-YDL 
Redstone  Arsenal,  AL  35808-5500 

1  Commander 

US  Army  Mobility  Equipment 
Research  &  Development 
Command 

ATTN:  DRDME-WC 
Fort  Belvoir,  VA  22060 

1  Commander 

US  Army  Mobility  Equipment  Cmd 
4300  Goodfellow  Blvd 
St.  Louis,  MO  63120 
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1  Commander 

US  Army  Tank  Automotive 
Command 

ATTN:  AMSTA-CG 
Warren,  MI  48090 

1  Commander 

US  Army  Tank  Automotive 
Command 

ATTN:  AMSTA-TSL 
Warren,  MI  48090 

1  US  Army  Armor  &  Engineer  Board 
ATTN:  STEBB-AD-S 
Fort  Knox,  KY  40121 

1  Project  Manager 
Tank  Systems 

ATTN:  AMCPM-GCM-SW 
Warren,  MI  48090 

1  Program  Manager 
Tank  Systems 
ATTN:  AMCPM-GCM-SA 
Warren,  MI  48090 

1  Project  Manager 

Fighting  Vehicle  Systems 
ATTN:  AMCPM-FVS 
Warren,  MI  48090 

1  Project  Manager 

Tank  Main  Armament  Systems 
ATTN:  AMCPM-TMA-PA 
Dover,  NJ  07801 
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Copies  Organization 

1  Commander 

Combined  Arms  Center 

ATTN:  ATZL-CG 

Fort  Leavenworth,  KS  66027 

1  Commander 

US  Army  Combined  Arms  Combat 
Developments  Activity 
ATTN:  ATCA-SA 
Fort  Leavenworth,  KS  66027 

1  Commander 

Concepts  Analysis  Agency 
8120  Woodmont  Avenue 
Bethesda,  MD  20014 

1  Commander 

US  Army  Foreign  Science  & 
Technology  Center 
ATTN:  AMXST-MC-3 
220  Seventh  Street,  NE 
Charlottesville,  VA  22901 

1  Commander 

US  Army  Training  &  Doctrine  Cmd 

ATTN:  ATCD 

Fort  Monroe,  VA  23351 

1  Commander 

US  Army  Training  &  Doctrine  Cmd 
ATTN:  ATCD-M 
Fort  Monroe,  VA  23351 
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1  Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SL,  Tech  Lib 
White  Sands  Missile  Range, 

NM  88002 

1  Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-TEM 
White  Sands  Missile  Range, 

NM  88002 

1  Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-TFB 
White  Sands  Missile  Range, 

NM  88002 

1  Commander 

US  Army  Field  Artillery  School 
Fort  Sill,  OK  73503 

1  Commandant 

US  Army  Infantry  School 
ATTN:  Infantry  Agency 
Fort  Benning,  GA  31005 

1  Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Benning,  GA  31005 


No.  of 

Copies  Organization 

1  Commander 

US  Army  Armor  Center 
ATTN:  ATZK-CG 
Fort  Knox,  KY  40121 

1  Commander 

US  Army  Armor  Center 
ATTN:  ATZK-CD 
Fort  Knox,  KY  40121 

1  Commander 

US  Army  Armor  Center 
ATTN:  ATZK-M1 
Fort  Knox,  KY  40121 

1  Commander 

US  Army  Development  & 
Employment 
Agency 

ATTN:  MODE-TED-SAB 
Fort  Lewis,  WA  08433 


1  Air  Force  Armament  Laboratory 
ATTN:  AFATL/DLODL 
Eglin  AFB,  FL  32542-5000 

1  Commander 
US  Marine  Corps 
ATTN:  APW 
Washington,  DC  20380 
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1  Commander 

Marine  Corps  Development  & 
Education  Command 
ATTN:  Plans  &  Studies  Div 
Quantico,  VA  22134 

1  Commander 

Naval  Weapons  Center 
China  Lake,  CA  03555 

1  Commander 

Marine  Corps  Development  & 
Education  Command 
ATTN:  Firepower  Division 
Quantico,  VA  22134 

1  AFELM,  The  Rand  Corporation 
ATTN:  Library-D 
1700  Main  Street 
Santa  Monica,  CA  00406 

Aberdeen  Proving  Ground 

Dir,  USAMSAA 

ATTN:  AMXSY-D 

AMXSY-MP  (H.  Cohen) 
Cdr,  USATECOM 

ATTN:  AMSTE-TO-F 
Cdr,  CRDC,  AMCCOM,  Bldg.  E3516 
ATTN:  SMCCR-RSP-A 
SMCCR-MU 
SMCCR-SPS-IL 
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